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ABSTRACT 



> ■ 

' Aims. In this paper we derive the structure of the Galactic stellar warp and flare. 
OQ . Methods. We use 2MASS red clump and red giant stars, selected at mean and fixed heliocentric distances of R e ^ 3, 7 and 17 kpc. 
CO ■ Results. Our results can be summarized as follows: (i) a clear stellar warp signature is derived for the 3 selected rings, proving that the warp 
C*~) ' starts already within the solar circle; (ii) the derived stellar warp is consistent (both in amplitude and phase-angle) with that for the Galactic 
interstellar dust and neutral atomic hydrogen; (iii) the consistency and regularity of the stellar-gaseous warp is traced out to about Roc ~ 20 
■ kpc; (iv) the Sun seems not to fall on the line of nodes. The stellar warp phase-angle orientation ((f> ~ 15°) is close to the orientation angle of 
the Galactic bar and this, most importantly, produces an asymmetric warp for the inner R a ^ 3 and 7 kpc rings; (v) a Northern/Southern warp 
C* ] , symmetry is observed only for the ring at R Q = 17 kpc, at which the dependency on i/> is weakened; (vi) treating a mixture of thin and thick 
O l' disk stellar populations, we trace the variation with Roc of the disk thickness (flaring) and derive an almost constant scale-height (~ 0.65 kpc) 
' within R G c ~ 15 kpc. Further out, the disk flaring increase gradually reaching a mean scale-height of ~ 1.5 kpc at R GC ~ 23 kpc; (vii) the 
> j . derived outer disk warping and flaring provide further robust evidence that there is no disk radial truncation at Roc ~ 14 kpc. 

Conclusions. In the particular case of the Canis Major (CMa) over-density we confirm its coincidence with the Southern stellar maximum warp 
occurring near / ~ 240° (for R B -7 kpc) which brings down the Milky Way mid-plane by ~ 3° in this direction. The regularity and consistency 
• • ■ of the stellar, gaseous and dust warp argues strongly against a recent merger scenario for Canis Major. We present evidence to conclude that 
all observed parameters (e.g. number density, radial velocities, proper motion etc) of CMa are consistent with it being a normal Milky Way 
outer-disk population, thereby leaving no justification for more complex interpretations of its origin. The present analysis or outer disk structure 
] ' does not provide a conclusive test of the structure or origin of the Monoceros Ring. Nevertheless, we show that a warped flared Milky Way 
contributes significantly at the locations of the Monoceros Ring. Comparison of outer Milky Way H i and CO properties with those of other 
galaxies favors the suggestion that complex structures close to planar in outer disks are common, and are a natural aspect of warped and flaring 
disks. 

Key words. Galaxy: structure - Galaxy: formation - galaxies: interactions - galaxies: dwarf - galaxies: individual Canis Major - galaxies: 
individual Monoceros Stream 



1. Introduction 

There was a revitalisation of interest in the outer disk of the 
Milky Way with the photometric discovery of significant over- 
densities of F stars in the Sloan Digital Sky Survey by Newberg 
et al. d2002> towards the Galactic anti-center. The reported 
over-densities defined a type of stellar Ring structure, rela- 
tively confined to the Galactic plane (\b\ < 30°), and stretching 
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* Based on archival data collected at the European Southern 
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over the sky for ~ 100° in longitude centered in Monoceros. 
A kinematic and spectroscopic study of the Ring structure by 
Yanny et al. ( t2003ll confirmed (i) a Galactocentric distance 
of 18 and 20 kpc (respectively above and below the Galactic 
plane); (ii) an inconsistency of its velocity dispersion with typ- 
ical Galactic halo and thick disk structures; and (iii) a metal- 
licity of [Fe/H]^ -1.6 consistent with halo populations. In the 
context of continuing searches for outer structure related to the 
Sgr dwarf, and other possible accretion events the most natu- 
ral scenario proposed by Yanny etal for the Monoceros Ring 
(Mon. Ring) was that it traced a remnant dwarf satellite galaxy 
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in the late process of disruption (if not already dissolved). 
Accumulating evidence supporting the existence of the Mon. 
Ring (Rocha-Pinto et al. l2l)0"3l Crane et al. 120031 Frinchaboy 
et al. [20041 soon triggered a search for a progenitor, obviously 
carried out near the Galactic plane. 

At the same time, in an alternative, more conservative sce- 
nario, the existence of the Mon. Ring was seen as the conse- 
quence of perturbations in the outer disk caused by ancient 
warps (Ibata et al. 2003}. Indeed, the Mon. Ring rotates in a 
prograde orbit that is almost circularized, strongly suggestive 
of a disk origin. Nevertheless, simulations by Helmi (|2003i 
showed that accretion models, where co-planar streams can 
follow circular orbits, were indeed feasible. This encouraged 
a search for a progenitor of the Ring: if there is one (the 
Sagittarius dwarf) why not another. 

In an analysis of 2MASS data, Martin et al. (I2004> as- 
sumed a symmetric Galactic vertical stellar distribution around 
b — 0° and searched for asymmetrical differences between 
Northern and Southern star-counts. Among other features, they 
pointed to an elliptical-shaped stellar over-density centered at 
(/, b) = (240°, -7°). They interpreted this over-density as the 
core of a satellite galaxy currently undergoing in-plane accre- 
tion, namely the Canis Major (CMa) dwarf spheroidal galaxy, 
the best Ring progenitor candidate. 

In Momany et al. d2004bt . we highlighted the fact that in 
the analyses in Martin et al. (20()4i and Bellazzini et al. (2004 
first astro-ph version) the possible influence of a Galactic stellar 
warp on the detection of a vertically asymmetric distribution so 
close to the Galactic plane was not considered. Observationally, 
the warp is a bending of the Galactic plane upwards in the first 
and second Galactic longitude quadrants (0° < / < 180°) and 
downward in the third and fourth quadrants (1 80° < / < 360°). 
We emphasise that the suggested reality of an outer stellar warp 
was not a new proposal by us. Among other earlier studies, one 
of particular relevance is that of Carney & Seitzer f 1993 1 who 
analysed the "Galaxy's own structure to obtain at least a peak 
at the outer disk", i.e. using certain lines of sights, one can 
look away from the plane, reducing the foreground disk sig- 
nal as well as reddening-absorption-crowding and derive the 
age and metallicity of the outer warped disk. Analyzing the 
color-magnitude diagrams of fields very near to the CMa cen- 
ter [(/, b) = (245°, -4°)] Carney & Seitzer claimed to have 
detected the main sequence and turnoff region of the outer 
Galactic disk. On the other hand, the analysis of Martin et 
al. (I2004> and Bellazzini et al. J2004> discounted the Galactic 
warp in this zone, so that the stellar populations previously 
identified as outer disk main sequence, were now proposed 
as an un-expected stellar population. Further hints of a star- 
count anomaly at I = 240° are found in Alard (120001 who, 
again, associated the "strong asymmetry" in this region with 
the Galactic stellar warp. 

In Momany et al. J2004b> we concluded that the CMa 
over-density can be fully accounted for if the Galactic disk (at 
/ = 240°) is 2° displaced/warped below the mid-plane; i.e. the 
symmetry axis for this region is at b — -2° and not b — 0°. In 
response to our analysis, Martin et al. (120041 presented radial 
velocities and argued that the Galactic stellar warp (location 
and amplitude), cannot explain the CMa over-density. In par- 



ticular, they argued that (i) a warp angle of -2° is not enough 
to erase the CMa over-density; and (ii) the CMa over-density 
is stronger in amplitude and located too far from the Southern 
hemisphere warp at / = 270°. More recently, Rocha-Pinto et al. 
(2005 1 proposed that the amplitude of the CMa over-density is 
small with respect to another over-density, this time located in 
Argo (I ~ 290°). In this later analysis, Rocha-Pinto et al. view 
the CMa over-density as the consequence of a dust extinction 
window aligned with a maximum warp location at I = 245° 
(as seen in Lopez-Corredoira et al. 120021 . Confusingly, Martin 
et al. ( 120041 and Bellazzini et al. ( 120061 use the same source 
(Lopez-Corredoira et al. 2002 1 to argue that the maximum warp 
location is at / = 270°. 

In order to clarify this situation, it is timely to re-evaluate 
the detailed properties of the Galactic stellar warp, specifically 
the location and amplitude of its maximum. Most recently, 
Conn et al. (I2005i presented a wide-field survey of the Mon. 
Ring and emphasised that the presence of Ring streams above 
and below the Galactic plane argue against a Galactic origin of 
the Mon. Ring. Moreover, they suggested that positive detec- 
tions of the Mon. Ring below the plane may also be correlated 
with the Triangulum- Andromeda (Tri And) structure (Majewski 
et al. 120041 and Rocha-Pinto et al. I2005i . In this scenario, al- 
though the Mon. Ring and TriAnd structure are located at dif- 
ferent distances, the TriAnd structure could be the distant arm 
of a multiply-wrapped tidal stream. 

In the main part of this paper we use 2MASS data to derive 
and trace the signature of the Galactic stellar warp as a func- 
tion of Galactic longitude. Our working hypothesis is simple: 
had the CMa over-density been due to an extra-Galactic accre- 
tion, one would expect it to appear as a distortion on top of a 
large-scale structure, that is the Galactic stellar warp. As we 
shall demonstrate, CMa as an over-density is easily accounted 
for as being the maximum Southern stellar warp. This, in our 
opinion, allows Occam's Razor to indicate a clear preference 
for a Galactic structural origin of the CMa over-density. We 
then discuss evidence which we believe associates the Mon. 
stream with the warped and flaring Galactic disk. In several 
Appendices we briefly review and comment on important as- 
pects of the observational properties of CMa and the Mon Ring. 

2. The Galactic Stellar Warp 

In Bellazzini et al. (2006i the authors commented on how the 
different adopted parametrization of the warp may lead to some 
confusion. Indeed, the warp has been derived from: (1) the 
mean latitude of the adopted tracer as a function of longi- 
tude (Djorgovski & Sosin .1989); (2) the ratio of star counts 
in Northern and Southern hemispheres as a function of lon- 
gitude (Lopez-Corredoira et al. [ZOO?*; and (3) the latitude of 
peak brightness as a function of longitude (Freudenreich et al. 
119941 . There has been some disagreement about these determi- 
nations. Bellazzini et al. ( 120061 argued that the Momany etal 
parametrization of the warp was "not a fair description" of the 
South/North over-densities. In the following we will firstly ex- 
pand on our method of tracing the Galactic warp and then sec- 
ondly we will compare our results with those derived by other 
methods. This comparison show excellent consistency. 
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In Momany et al. J2004bl Fig. 2) we extracted dereddened 
2MASS M-giants (0.85 < (J-K) < 1.3) within an oblique box 
surrounding the CMa red giant branch in the CMD, and falling 
in a strip between -20° < b < 20° and 235° < / < 245°. 
Initially assuming reflection symmetry around b — 0° (warp an- 
gle w - 0°) we performed Northern and Southern star counts, 
binned in 0.5°. This showed clearly (upper panels of the same 
Figure) that by assuming a symmetry around b = 0°, as done in 
Martin et al. (2004). one can recover the identified (CMa) over- 
density. Successively the symmetry axis was varied within the 
range \b\ < 5°, in steps of 0.1°. For each step (i.e. each warp 
angle) the Northern and Southern latitude profiles were folded 
and a reduced^ 2 defined as: 

X 2 ~ / \[(y North - y South) 1 1 iy North + ySouth)] (1) 

was employed to determine the warp angle which minimized 
the differences between the two profiles (y is the Log(N) of star 
counts computed in the latitude profiles, see Fig.0}. A symme- 
try of around b = -2° (i.e warp angle w = -2°) almost com- 
pletely canceled the (CMa) over-density. This warp angle has 
been derived [not modeled as argued in Martin et al. (2004) and 
Bellazzini et al. ( I2004H and agrees with the angle derived for 
the gas warp by Freudenreich ( 1994). It is important to note that 
this result is obtained independently of the adopted deredden- 
ing method: (i) pure Schlegel et al. J1998b values; or (ii) modi- 
fied with the Bonifacio et al. (2()()0i correction [see Amores & 
Lepine i2QQ5| for recent confirmation of errors in the Schlegel 
et al. values]. 

Thus, we are able to measure the latitude angle at which the 
latitude profiles crossing the Galactic disk are symmetric; i.e. 
the warp-angle or the mid-plane of the warped disk. In this sec- 
tion we apply this method for the entire 2MASS catalog within 
-20° < b < 20°, and derive the global Galactic warp as a 
function of longitude. Most importantly, the warp signature is 
traced by means of different stellar sources [red clump (RC) 
and red giants (RGB)] at different distances. This is particu- 
larly important because we want to investigate the impact of 
the tracer contamination and its distance on the derived warp 
signature. 

To emphasize the particular importance of the contamina- 
tion and distance properties of the assumed tracer (and the in- 
trinsic difficulty for similar investigations) we will briefly com- 
ment on two recent studies probing the Galactic warp, namely 
Yusifov (2004) and Lopez-Corredoira et al. ( I2002> . 

2.1. The Yusifov stellar warp model 

In his analysis, Yusifov used the asymmetric distribution of 
1412 Galactic pulsars (from the Manchester et al. 120051 cata- 
log) to trace the Galactic stellar warp and flare (the increase in 
scale-height as a function of Galactocentric distance). Yusifov 
limited his analysis to pulsars within \b\ < 5° whose distances 
are > 1 kpc. He calculated the ratio of the cumulative number 
of pulsars above and below the Galactic plane as a function of 
Galactic longitude, and so derived a warp model. However, as 
seen in Fig.^ the Manchester et al. ( 120051 ) catalog has only 1 
pulsar at R > 5 kpc between 200° < / < 270°. Clearly, the pul- 
sar catalog is already strongly incomplete at R ~ 5 kpc in the 
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Fig. 1. The distribution pulsars in the Manchester et al. (|2005! 
catalog. Filled squares and open circles are pulsars with 
Galactocentric distances between 1 and 5 kpc, and 5 and 10 
kpc respectively. Open triangles are pulsars at Rqc ^ 10 kpc. 

outer Milky Way. Although one might doubt how these data 
can reliably predict the number density of stars at CMa dis- 
tances, it is true that the Yusifov warp model predictions have 
turned out to be comparable with other studies based on more 
complete samples, and we therefore will compare our results 
with this model. 

2.2. The Lopez-Corredoira et al. stellar warp model 

The Lopez-Corredoira et al. 2MASS-based investigation was a 
major attempt at studying the Galactic stellar warp and flare. 
There are, however, a few points that must be kept in mind re- 
garding their analysis. Firstly, 2MASS was not complete at the 
time that their analysis was concluded, so that they were able 
to analyse only 820 lines of sight, each of area between 0.5 
and 1.0 square degrees, restricted to \b\ = +3°, +6°, +9°. This 
is particularly important in understanding the determination (cf 
Martin et al. l^UOll and Bellazzini et al.E001|2QQ6| of the maxi- 
mum of the stellar warp at / ~ 270°, a result often cited from the 
Lopez-Corredoira et al. paper. Indeed, the area around I ~ 270° 
was missing in the Lopez-Corredoira et al. study, and nowhere 
in that article is it stated that the stellar warp maximum is at 
I ~ 270°. Indeed the formula describing the Lopez-Corredoira 
et al. warp model at CMa distances shows that the maximum 
stellar warp is near to I ~ 240° (see Section 0}. The value of 
I ~ 240° for the maximum Galactic warp is in fact found in 
other studies, e.g. Freudenreich ( 119941 their Fig. 3) studying 
the gaseous warp, and Djorgovski & Sosin ( 119891 their Fig. 1) 
analyzing IRAS sources 1 . 

Second, Lopez-Corredoira et al. derive the ratio of Northern 
to Southern star-counts, Rrc = NNorthlNsouth for red clump 
stars with K < 14.0 (see their Fig. 15). To help us evaluate 
the impact of contamination when using RC stars as warp- 
tracers, in Fig. 13 we use Besangon (Robin et al. I2003i simu- 
lated color-magnitude diagrams around CMa and plot the data 
as a function of their distance. Focusing our attention around 
(J - K) ^ 0.6, one sees that the stellar populations at < 2 kpc 

1 In particular note a clear warp excess around / ~ 240° with respect 
to their fitted function. 
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Fig. 2. The upper 4 panels display a Besancon 1° square 
area simulated color-magnitude diagram around CMa (I, b)- 
(240°, -7°). These show the expected different contributions 
of stars as a function of their distance. The lowest panel dis- 
plays data from 2MASS between 238° < / < 242° and 
-8° < b < -6°. The outlined vertical boxes mark 4 distance 
intervals as selected from the red clump distribution in simu- 
lated diagrams. The 3 oblique boxes trace 3 samples of red gi- 
ant stars extracted at three different heliocentric distances. The 
thick line is the Majewski et al. calibration of the Sagittarius 
dwarf RGB, shifted to 2.8, 7.3 and 16.6 kpc. 




Kiloparsec 



Fig. 3. A schematic view (from Drimmel & Spergel 20QD of 
the Milky Way as seen from its North pole showing the 4 spi- 
ral arms as mapped by H n regions and dust (see also Vallee 
2005 and Russeil 2003). The Galactic center, the Sun and the 
names of the spiral arms are plotted. Note the presence of the 
Local arm (Orion-Cygnus) close to the Sun's position and the 
outer arm that is also called Norma-Cygnus. The 3 heliocentric 
circles define the regions at fixed R Q where we extract RGB 
samples to trace the Galactic warp and flare. The direction of 
the warp phase angle <p, and longitude directions are also plot- 
ted. We have also added the Galactic bar at an orientation an- 
gle of 14° (following Freudenreich [1998 1 with length 3 kpc 
(following Vallee 2005 1. Note however that a recent Galactic 
Legacy Mid-Plane Survey Extraordinaire (GLIMPSE) study by 



Benjamin et al. 1 2005 ) report on a linear Galactic bar with half- 
length of4.4 + 0.5kpc. 



which will contribute to the final R^f- ratio include: (i) RC 
stars at 8.0 < K < 10.0, and most importantly, (ii) dwarfs 
at the fainter magnitudes K a > 11.5. At magnitudes around 
the RC of CMa (K ~ 13.0), we estimate a contamination by 
local foreground dwarfs of about 25-30%. Thus, if one is inter- 
ested in the Rrc ratio at say, the CMa distance, one must bear 
in mind that this ratio is subject to potentially serious contami- 
nation by RC and dwarf stars at distances closer than the CMa 
over-density. 

2.3. Tracing the stellar warp and flare with RC and 
RGB stars 

We employ our method in tracing the Galactic warp and flare 
using two stellar tracers: red clump and red giant stars. For both 
tracers we repeat the analysis for different mean distances. The 
advantage in using two stellar tracers derives from their dif- 
ferent contamination status. With respect to RC stars, red gi- 
ants are an ideal instrument in probing the Galactic warp and 



flare since these suffer less external contamination by nearby 
dwarfs. At the same time red giants projected at different dis- 
tances do not overlap in the color-magnitude diagram (i.e. no 
internal contamination) and this guarantees a better distance 
separation. Indeed, being the bright and evolved part of the 
faint and un-detected main sequence stars, red giants allow us 
(already in the not so deep 2MASS catalog) to probe the most 
distant, and almost entire, Milky Way disk populations. 

The vertical boxes in Fig. |2 set 4 heliocentric distance in- 
tervals (< 2, 2 -r 4, 4 -r 6 and 6 4- 8 kpc) from which we extract 
RC stars and use their star-counts to derive the warp and flare. 

Similarly, the 3 oblique boxes trace red giant stars at 
fixed heliocentric distances between 2.3 -f 3.2, 5.5 4- 9.1, and 
13.2 4- 20.0 kpc. Thus, the mean assumed distances of the 
RGB samples are 2.8, 7.3 and 16.6 kpc, with the intermedi- 
ate sample being centered on the CMa RGB. To estimate the 
distances to Milky Way M-giants we apply the same method 
as in Martin et al. (2004). that first led to the identification of 
the CMa over-density. This is done by using the Sagittarius 
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Fig. 4. The upper panels display two examples of the deriva- 
tion of best warp angles and their relative error by means of 
a reduced x 1 - The lower panels show the latitude profiles for 
a third field: (i) assuming no warp (upper panel) and; (ii) hav- 
ing derived the best warp angle (lower panel). The dashed and 
solid lines refer to the Northern and Southern latitude profiles, 
respectively. All panels refer to the RGB sample at R Q = 7.3 
kpc. 

dwarf RGB as a reference. Assuming a distance modulus of 
16.9 and a mean [Fe/H] of ~ -0.5 2 Majewski et al. (120031 de- 
rived the following calibration of the Sagittarius RGB: K = 
-8.650 x (J - K) + 20.374 (see the thick line in Fig.EJ. The 
mean metallicity of Sagittarius can be considered intermedi- 
ate between inner (more metal-rich) and outer (more metal- 
poor) Milky Way disk stars. Indeed, the mean abundance of the 
disk stars is known to vary between -1.0 <[Fe/H]< +0.3 [see 
Bensby, Feltzing & Lundstrom d2004> for differences between 
thin and thick disk populations] 

2 A recent study by Monaco et al. I2005i report on high resolution 
UVES spectroscopy showing a mean metallicity of [Fe/H]= -0.41 + 
0.20. 



Fig. 5. An example of deriving the scale height (Zj,) from the 
RGB sample at R Q = 16.6 kpc. The filled circles show the 
Northern hemisphere latitude profile while the open circles 
show the Southern profile. Vertical lines delimit the region 
where we apply a linear fit (dotted line) to the data. 



At the same time we note that the Majewski et al. RGB cali- 
bration would reproduce the entire range of disk metallicities if 
age is allowed to vary. For example the Sagittarius RGB would 
overlap with a 10 Gyr and [Fe/H]= -0.7 theoretical isochrone 
(appropriate for outer disk populations?) as well as one with 
4 Gyr and [Fe/H]= -0.4 (appropriate for inner disk popula- 
tions?). A systematic uncertainty in estimating the distances 
of the Galactic M-giants is therefore un-avoidable. Similarly, 
contamination between different populations and uncertainties 
in reddening corrections will cloud our analysis. A schematic 
view of the Milky Way and of the regions where we probe the 
stellar warp and flare are shown in Fig. [3] 

Examples of the application of our method are shown in the 
upper horizontal panels of Fig.|4] reporting the derivation of the 
best warp angle for two fields extracted from the R Q = 7.3 kpc 
RGB sample. We show two cases: the first field 80° < / < 90° 
and \b\ < 20° is more populated and the warp angle (minimum 
X 2 , as reported in equation^ is easier to determine, having a 
relatively small 0.2° error. The second field, 240° < / < 250°, 
is less populated and the error is higher reaching 0.5°. 

The lower panels of Fig.@J show an example of deriving the 
best warp angle for latitude profiles between 220° < / < 230° 
from the R = 7.3 kpc RGB sample. The upper panel shows 
the Northern and Southern latitude profile assuming no warp, 
i.e. a symmetry axis around b = 0°. The lower panel shows the 
same two profiles assuming a symmetry axis around b = -2.9°, 
as derived when the differences between the folded profiles are 
minimized and the symmetry angle is allowed to vary within 
\b\ < 5°. One may note that for large warp angles the Southern 
latitude profile becomes shorter than the Northern one. This is 
due to the warp angle becoming significant compared to the 
imposed angular limit in the extraction; i.e. \b\ < 20°. 

Once the warp angle has been derived for each line of sight 
we proceed in estimating the vertical density profile of the 
stellar disk: the scale-height (Z/,)- Our intent is to character- 
ize the radial trend of the stellar disk scale-height. A disk flar- 
ing is usually seen as an increasing scale-height towards the 
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outer parts of the disk. An example of the adopted procedure is 
shown in Figure[5] Having derived the warp angle, the b angles 
were transformed in linear Z height according to the distance 
from the Sun of the sample. Then the Northern and Southern 
profiles were overlapped and the star density profile was fitted 
by a power law, thus deriving Z/,. The fitting of the vertical den- 
sity profile was made within 2 limits in Z, so as to exclude: (i) 
the highly obscured inner regions within +2.0° from the mid- 
plane and; (ii) the very outer regions where in cases of high 
warp angles the latitude profile of one hemisphere is short with 
respect to the other hemisphere (cf the latitude profiles in lower 
panel of Fig.|4}. 

The vertical Z scale-height was fitted by a single exponen- 
tial. We have not tried to use a more complicated formula like 
the sech or double exponential (see for example the analysis 
by Alard 200(3. Thus, our analysis is aimed at measuring the 
order of magnitude of the flare in the outer disk. A detailed 
parametrization of the flare (in terms of thin-thick disk sepa- 
ration), although very interesting, is beyond the scope of this 
paper. In particular, we note that for the regions of most inter- 
est (the outer disk) any thin disk flaring would act in a way that 
increases the confusion between the thin and thick disk compo- 
nents. This in-ability in separating thin-thick disk components 
in the outer regions is exacerbated by the unavailability of kine- 
matic all-sky data, which might allow a separation of the two 
components. Thus, for the outer disk, one is left with determin- 
ing only an approximate amplitude of any flare. Figure[5]shows 
an example of the vertical density profile fit for the RGB sam- 
ple at a heliocentric distance of 16.6 kpc (Rcc - 14.9 kpc). The 
abscissa, |Z - w\, indicates the linear Z height folded around the 
mid-plane of the disk as determined by the warp angle analysis. 
The analysis of the radial trend of the vertical scale height, Z/,, 
for the RGB sample will be discussed in Section 5. 

2.4. Comparing DIRBE integrated surface photometry 
with 2MASS star counts 

In this paper we compare the derived stellar warp with that 
obtained for other Galactic components, including integrated 
light, neutral gas and the interstellar dust. In regards to this, we 



mainly make use of the Freudenreich et al. (1994) study where 
the warp due to dust and integrated light has been derived using 
the Diffuse Infrared Background Experiment (DIRBE) map- 
ping of the Galactic plane (see also Vig, Ghosh & Oiha l2005Y 

The results of DIRBE however come from integrated light 
surface photometry which sum contributions from composite 
stellar populations at different distances and luminosities, and 
therefore these can be correctly compared with the 2MASS 
results only after a discussion of their relative weight. The 
DIRBE 1.2/im and 2.2/im spectral bands correspond roughly 
to the infrared J and K pass-bands. It is important to note that 
the DIRBE near-infrared (A < 5//m) emission is dominated by 
stellar disk stars [whereas in far-infrared bands (e.g. 240jt/m) 
the emission by interstellar dust dominates]. 

To better understand which stellar populations contribute 
most in the DIRBE emission maps, we used the 2MASS star 
counts to integrate the light coming from specific stellar pop- 



ulations. We limit this analysis to only one line of sight (the 
2MASS color-magnitude diagram presented in the lowest panel 
of Fig. |3, and expect this to be representative of any other 
line of sight. We roughly disentangle three main populations 
in this 2MASS color magnitude diagram: (i) MS objects with 
(J - K) < 0.45, mainly main sequence and blue super-giants 
stars; (ii) RC stars with 0.45 < (J - K) < 0.65; and (iii) RGB, 
red super-giant and asymptotic branch stars with (J-K) > 0.65. 
For each of these three stellar populations we summed up the 
light contribution coming from the single point sources (stars) 
falling within the above color intervals. 
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Fig. 6. The relative luminosity contribution of the three stellar 
populations in the J band (upper panel) and K band (middle 
panel). The lower Figure shows the cumulative distributions of 
the distances of RGB stars without any luminosity-weighting 
along the line of sight (dark line). The grey line shows instead 
the distance-weighted distribution. 



In Fig. [6] we show the relative luminosity contribution of 
the three stellar populations in the J band (upper panel) and 
K band (middle panel). A limiting magnitude of K — 14.0 has 
been applied in estimating the cumulative and relative luminos- 
ity of the three different stellar populations. The middle panel, 
displaying the relative luminosity contribution of the three pop- 
ulations to the total luminosity, shows clearly that in the case 
of the K band the total light contribution of MS stars is almost 
negligible in DIRBE, contributing only 15% of the total light. 
On the other hand, RGB stars contribute up to 60% of the total 
light, a level that goes up to 82% once also the RC are summed 
together with RGB stars. Thus we conclude that the level of 
contamination by MS stars in 2.2yum (~K) DIRBE emission 
maps is fairly low and that red populations dominate the inte- 
grated light. Similar conclusions can be drawn for the shortest 
DIRBE wavelength band at 1 .2//m (~J, upper panel): the RGB 
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and RC populations dominate the DIRBE light, contributing up 
to 75% of the total light. 

We also checked the weighted distance range at which 
DIRBE is more sensitive. Considering that the RGB population 
is the dominant population in both J and K bands we calculated 
their distance distribution using the Majewski calibration of the 
Sagittarius red giant branch (see Sect. 12. 3> . The dark line in the 
lowest plot of Fig.[6]shows the cumulative distributions of the 
distances of RGB stars without any luminosity weighting along 
the line of sight. The grey line shows the cumulative distance 
distribution once we weight the stellar distance with the stars 
light. 

Clearly, the DIRBE integrated light is linearly sensitive to 
distances within ~ 9 kpc (at least in this direction I ~ 240°) 
where the light contribution reaches 80% of the total. This re- 
sult, based on RGB stars, should not change considerably once 
we add the contribution of (i) the RC stars which weight more 
to larger distances; and (ii) the MS stars which weight more 
to shorter distances. Thus, being of the same size, the two op- 
posite contributions of the MS and RC stars will to first order 
cancel out. 

In conclusion, with the help of the 2MASS star counts 
we have identified which stellar populations are the main con- 
tributers to the DIRBE luminosity, and at what distances these 
contribute the most. The results are that the RGB stars are the 
dominant contributors of the J and K DIRBE integrated light, 
within ~ 9 kpc. This allows us to perform a fruitful compar- 
ison of the DIRBE warp results with those obtained from our 
2MASS analysis. 
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Fig. 7. The stellar warp as derived from RC stars at R Q < 2 kpc 
(filled squares). The results are compared with the latitude of 
peak brightness obtained from DIRBE at wavelengths 1.2, 3.5 
and 240 pm by Freudenreich et al. J 19941 grey open triangles). 



3. The stellar warp as traced by RC stars 
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Fig. 8. The stellar warp as derived from RC stars for 4 distance 
intervals (dark filled squares). The results are compared with 
the latitude of peak brightness obtained from DIRBE at wave- 
lengths 3.5 and 240 pm by Freudenreich et al. J 19941 grey open 
triangles). 

In this section we present the derived stellar warp as traced 
by RC stars. We anticipate two obvious points: (i) the warp 
amplitude depends on the employed wavelength; shorter wave- 
length investigations probe nearer areas and therefore may un- 
derestimate the warp amplitude if it increases with distance; 
and (ii) the warp amplitude derived depends on the contamina- 
tion of the employed stellar tracer: inclusion of nearby stars in 
samples of more distant stars will underestimate the warp am- 
plitude at larger distances. Overall, red clump stars are not the 
ideal tracers of the warp at larger distances, and we will show 
this using Fig.0and [H] 

The upper left panel of Fig. displays the stellar warp 
as derived using RC stars at R < 2 kpc. At first sight, the 
presence of a global stellar warp can be debatable, as it con- 
trasts with the neat sinusoidal function seen in Djorgovski & 
Sosin (1989 1 for IRAS sources. It remains true however that the 
Southern stellar warp is distinguishable between 190° -r 290°. 
This brings the Galactic mid-plane ~ 2° down at around / ~ 
250°, and demonstrates that the stellar warp is detectable al- 
ready in the solar circle. Problems arise when searching for a 
global stellar warp in the Northern hemisphere. A sudden drop 
at around I ~ 120° seems to interrupt a global warp signa- 
ture. To understand the absence of a clear stellar warp in the 
Northern hemisphere we make a comparison with the latitude 
of peak brightness as derived in Freudenreich et al. I1994t us- 
ing DIRBE mapping of the Galactic plane. Overall, the DIRBE 
data-points at 1 . 2pm reproduce the global features derived from 
nearby R < 2 kpc RC stars. In particular both data-sets 
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show an early drop at / ~ 120° demonstrating: (i) the consis- 
tency of the DIRBE-2MASS comparison; and (ii) the presence 
of a localized nearby structure. Indeed, the drop in 2MASS 
and DIRBE traces the presence of the Orion-Cygnus segment 
(called Local Arm in Russeil l2003> . 

The lower left panel of Fig. over-plots the DIRBE data- 
points at 3.5pm. As anticipated before, at A < 5fim DIRBE 
emission is still dominated by disk stars. However, as one 
moves to longer wavelengths DIRBE data-points sample more 
distant structures that were obscured at shorter wavelengths. 
This explains why at 3.5/im we do not observe the strong drop 
at / ~ 120° (seen at l.2pm and in RC stars at R e < 2 kpc) 
and this is replaced by a global, and smooth warp signature. 
This effect is further demonstrated when comparing the warp 
obtained for RC stars at R s < 2 kpc with that due to the 
dust (traced at 240jum). The dust warp does not show the drop 
at I ~ 120° anymore, and although still fluctuating, overall a 
global and large-scale structure (the warp) is identified. We em- 
phasize however that a simple comparison between the warp of 
RC stars at R G < 2 kpc with that at 240yum is rather improper, 
since the DIRBE dust warp is sensitive to more distant regions. 
The comparison has been meant to show how increasing the 
adopted wavelength unveils more distant regions and a larger 
warp amplitude. 

Before drawing our conclusions on the RC stars as a stellar 
warp tracer, we analyze the impact of their contamination sta- 
tus on the global warp properties. This is illustrated in Fig. [8] 
which displays the stellar warp as derived from RC stars at he- 
liocentric distances between 2+4,4 + 6,6 + 8 and < 8 kpc. The 
upper panels of Fig. Improve again how the presence of a global 
warp can be debatable for the "nearest" RC samples. Things 
change considerably for RC stars between 6 < R Q < 8. At this 
distance a global warp signature is evident in both hemispheres. 
The warp is regular and shows a strong asymmetry in the maxi- 
mum warp amplitude. At I 210° the warp brings the Galactic 
mid-plane ~ 2.7° below the nominal b = 0°, whereas in the 
Northern hemisphere the maximum warp amplitude stops at 
b = +1.5°. Comparing the warp obtained for the 6 < R B < 8 
kpc RC sample with DIRBE data shows excellent agreement 
with that obtained at 240yum tracing the dust. This indicates 
that the RC star sample between 6 < R Q < 8 kpc (although it 
includes up to 30% contamination by local dwarfs) still enables 
sufficient isolation of distant stars, those whose warp signature 
matches that of DIRBE at 240jum. 

The impact of contamination on the retrieved warp ampli- 
tude is best illustrated in the lower right panel of Fig. [H] With 
respect to the warp derived from RC stars between 6 < R e < 8 
kpc, those between < R Q < 8 kpc show a significant change 
in the amplitude of the warp maximum and location. In par- 
ticular, the maximum warp amplitude passes from ~ -2.7° 
(for RC between 6 < R B < 8 kpc) to ~ -1.5° (for RC at 
< R < 8 kpc). A comparison with DIRBE data now shows 
a better agreement with the 3.5yum data and this is again the 
expected result. Indeed, both the 2MASS < R e < 8 kpc star 
counts and the DIRBE 3.5yum surface photometry are sensitive 
all along the line of sight up to ~ 8 kpc. 

In conclusion, the plots presented in Figs. Q 00 highlight 
the fact that a thorough characterization of the Galactic stellar 



warp, in particular its maximum amplitude and maximum lo- 
cation, sensitively depend on the distance and contamination- 
status of the employed tracer. The lower panels of Fig. [8] are 
the most significant in demonstrating how the maximum stel- 
lar warp can shift its location by ~ 40° in longitude and can 
halve its amplitude depending on the contamination degree of 
the adopted tracer. In addition, the lower right panel of Fig. [8] 
proves the validity of our method in tracing the stellar warp, 
and the consistency of the comparison between DIRBE and 
2MASS. 

This being established we will now focus our discussion on 
the stellar warp as derived using red giants. As anticipated in 
Sect. 12.31 the red giant samples include less contamination by 
local dwarfs and enable a better and more reliable separation of 
more distant stars. 

4. The stellar warp as traced by RGB stars 

Figures[9]and ^|display the stellar warp as derived from RGB 
samples at R e = 7.3 and R Q = 16.6 kpc. A global regular warp 
signature is clearly evident in both samples, reflecting large- 
scale Galactic structure (see Section l4~2t . The most striking 
feature for the RGB Rq = 7.3 kpc sample is the amplitude of 
the Southern maximum warp, which lowers the Galactic mid- 
plane by almost 3°. By applying a Gaussian fit to the data- 
points between 160° < I < 310°, we estimate a mean warp 
maximum at I ~ 235°. Similarly, for the R Q = 16.6 kpc RGB 
sample the warp maximum is estimated at I ~ 240°. 

The lower left panel of Fig. [9] shows a comparison be- 
tween the stellar warp as derived from the RC star samples 
(6 < R @ < 8 and < R Q < 8 kpc) with the RGB sample 
at R e = 7.3 kpc. This comparison is an excellent demonstra- 
tion of how (i) the distance and (ii) the contamination status of 
the adopted stellar tracer affect the discussion on the location 
and amplitude of the maximum stellar warp. The RC sample 
between 6 < R Q < 8 kpc shows a comparable Southern warp 
maximum amplitude with the RGB sample at R Q = 7.3 kpc, 
and this reflects the fact that both samples refer to the same dis- 
tance range. Yet, the location of the warp maximum for the RC 
sample between 6 < R < 8 kpc is ~ 20° shifted with respect 
to that for the RGB sample at R e = 7.3 kpc. On the other hand, 
the more contaminated RC sample between < R Q < 8 kpc 
shows a better agreement with the RGB sample at R Q = 7.3 kpc 
in the mean location of the warp maximum around / ~ 240°. 
However, the maximum amplitude of the RC sample between 
< R e < 8 kpc has clearly decreased, a reminder of how 
the inclusion of nearby stars leads to an under-estimated warp 
amplitude. 

We now turn our attention to the comparison of the stellar 
warp (derived from the RGB sample at R Q = 7.3 kpc) with 
that obtained for the interstellar dust and neutral atomic hydro- 
gen components. The upper right panel shows the comparison 
with the DIRBE data at 240/im. The agreement between the 
stellar and dust warp is evident. The two share the same warp 
phase-angle and show only small differences in the Southern 
maximum amplitude. 

A recent upgrade of the Drimmel & Spergel (1200 It 
Galactic dust model, using DIRBE 240/L/m data, shows the 
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Fig. 9. The upper left panel displays the stellar warp as derived from red giants at R Q 7.3 kpc. In the lower left panel we 
compare the derived warp with those based on two RC samples between: (i) 6 < R Q < 8, and (ii) < R Q < 8 kpc. The upper 
right panel over-plots the dust warp as derived from DIRBE at 240/mi data, whereas the lower right panel over-plots the Galactic 
warp as derived from neutral Hi gas (also from Freudenreich et al. 1 19941 1, The lower panel shows the density maps and contours 
of the R Q 7.3 RGB sample. One can trace by eye the warp signature (a colored version is more appropriate). Over-plotted is 
also the location of the mean mid-plane warped disk as a function of longitude. 
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presence of an extended spiral arm entering the third quad- 
rant and intersecting the suggested location of CMa. The re- 
construction of the Galactic dust distribution is always model 
dependent. However, the new model (Drimmel 2005 1 indicates 
that the far infrared DIRBE 240/im data, are sensitive to at least 
~ 7 kpc from the Sun, and possibly beyond. The distance con- 
sistency of the far infrared DIRBE 240/im data with the RGB 
star sample at R = 7.3 kpc (and given the excellent agreement 
presented in the upper right panel of Fig. |5) allows us to as- 
certain that the dust and stars at CMa distances are similarly 
warped. 

Besides the excellent agreement between the dust and stel- 
lar warp, a similar conclusion is also obtained for the gaseous 
and stellar warp. The H i latitude of peak brightness as derived 
by Freudenreich et al. ( }1994f ) follows the same warp signature 
we have obtained from the RGB sample at CMa distance. The 
only obvious difference is related to the maximum Southern 
warp amplitude of H i that seems to saturate; showing an al- 
most constant w warp-angle of -1.75° for longitudes between 
210° and 270°. Might this observed difference between the am- 
plitude of the stellar and gaseous warp suggest an "accretion- 
perturbation" on top of the large-scale warp structure derived in 
this paper?. In regard to the reality of this effect, one must keep 
in mind that over the last few years there has been alternating 
evidence on warp amplitude differences between stars and gas. 
Indeed, we remind the reader that the disk radial truncation 
hypothesis at Rqc — 14 kpc was first proposed (Freudenreich 
et al. 119941 and Porcel, Battaner & Jimenez- Vicente W 997\ in 
order to explain "why the stellar warp seemed half the ampli- 
tude of the gaseous warp" . Thus, our results which go in the 
opposite direction (stellar RGB warp at R Q 7.3 kpc showing 
twice the H i warp amplitude) should only be read in terms of 
the different confusion status, and therefore of the mean probed 
distances, of the gas. 

Indeed, whereas we are able to isolate the warp signature 
due only to the stellar populations at R G ^7.3 kpc, the gaseous 
warp derived by Freudenreich et al. i 19941 is that due to the 
summing of different warp signatures by gas distributed along 
the line of sight, at different distances. This is a direct con- 
sequence of the inability to recover gas distribution distances 
in lines of sight without clear kinematical signatures. Thus the 
saturation effect seen in the H i warp may reflect that seen for 
the stellar warp (cf the lower panels in Fig. [8] and the lower 
left panel in Fig. [9} when including nearby stars in distant star 
samples. More fundamentally, leaving apart the flattening of 
the gas warp in the range 210° < I < 250°, one notes excellent 
agreement between of the gas and stellar warp profiles at all 
the other viewing angles. 

Moving to the stellar warp derived by RGB stars at R Q ss 
16.6 kpc (Fig. 1 1 21 one notes structure ("swelling") between 
140° < I < 180°. Is this a deviation from a smooth global 
warp signature?. We suggest that this is the same effect seen 
for nearby RC samples, and it is due to the Orion-Cygnus Local 

3 Note that although the Freudenreich et al. paper was based on 
DIRBE data, they also derived the gaseous warp based on a re- 
analysis of the velocity-integrated HI emission maps from the survey 
of Weaver & Williams 115751 . 
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Fig. 10. Same as Fig |9] but for RGB stars at R G 16.6 kpc. 
Note the structure ("swelling") around I ~ 160° and how the 
the agreement with Hi is better with respect to the R Q 7.3 
kpc sample. 



Arm (short dashed line in Fig. [3j. This "swelling" is not seen 
in the R Q =a 7.3 kpc sample because the selected RGB stars 
are bright enough to avoid the faintest regions of the diagrams 
(mostly populated by nearby dwarfs). On the contrary, the RGB 
sample at R 16.6 kpc can be easily contaminated by nearby 
R < 2 kpc dwarf stars [compare the upper and lower panels of 
Fig. |2) which would enter our selection box if subject to high 
reddening or photometric error. In turn, the structure appears in 
the R Q 2.8 kpc RGB sample (Fig.ll2> supporting our inter- 
pretation that the "swelling" is due to nearby dwarf stars that 
contaminate primarily the R Q =a 16.6 kpc RGB sample. Given 
the similarities between the R Q 2.8 and 16.6 kpc warps, to 
what extent might the R e 16.6 kpc warp signature be in fact 
a signature of nearby contamination in the R Q 16.6 kpc RGB 
sample? The amplitude difference in the Northern warp of the 
R e 2.8 and 16.6 kpc samples however argues against this 
possibility (compare Figs.ll2land II 31 . 

Recently McClure-Griffiths et al. (I2004> presented an H i 
study from the Southern Galactic Plane Survey (SGPS), and 
pointed to the possible presence of a distant spiral arm in 
the fourth quadrant of the Milky Way. The distinct and cohe- 
sive feature has been traced over 70° and is located between 
18 < Rqc ^ 24 kpc. In regards to this, it is interesting to 
note that our detection of a stellar warp in the RGB sample at 
R - 16.6 kpc might represent the stellar counter-part of this 
distant H i spiral arm. Indeed, the RGB sample at a mean dis- 
tance of R G =! 16.6 kpc probes the stellar populations between 
13.2 -r 20.0 kpc from the Sun (Fig. [3} and this is compatible 
with the mean position of this H i arm being within ~ 20 kpc 
from the Galactic center (cf Figure. 3a of McClure-Griffiths et 
al. 120041 . 
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Interestingly, McClure-Griffiths et al. find that this distant 
arm is well-confined to the Galactic plane, dropping at most 
by ~ 1 kpc below the Galactic equator. This again is in agree- 
ment with our result (Fig. 1 1 31 that the Galactic mid-plane (at 
these distances and lines of sight) is located at ~ 0.5 kpc be- 
low latitudes of b — 0°. The SGPS H i survey analysis however 
remains limited to / ~ 250° (McClure-Griffiths et al. 120051 . 
so that a detailed comparison of the stellar-gaseous warp and 
the interpretation of this part of the outer disk awaits the com- 
pletion of the third quadrant. In conclusion, and as we shall 
argue in Section [5] the detection of the warp in the R Q 16.6 
kpc sample proves that the MW stellar disk is not truncated at 
Rgc — 14 kpc, and that a more extended stellar component is 
present. 

4.1. Comparison with warp models 

Before discussing the possible implications of the warp as de- 
termined empirically here on the CMa over-density, we first 
compare our results with other available models. In Fig. II llwe 
show a comparison of our derived Z - warp at the proposed 
distance of CMa (R Q ^7.3 kpc) with available models based 
on stellar tracers, namely Drimmel & Spergel GOOD 4 . Lopez - 
Corredoira et al. (12002b and Yusifov d2004> . The 3 models used 
the distribution of integrated star light, red clump and pulsar 
stars, respectively. In the following all models were converted 
to R Q = 8.5 kpc. At first sight, one sees how all three models 
converge on showing a warp maximum near I = 240° rather 
than / = 270°, a clear hint of the association of CMa with the 
warp maximum amplitude. Most importantly, at a heliocentric 
distance of 7.3 kpc the 3 models show excellent agreement with 
our observed warp amplitude around I ~ 240°. 

One also notes that there exists: (i) an overall good agree- 
ment between the Lopez-Corredoira et al. ( 2002 1 and Drimmel 
& Spergel (1200 1> models; (ii) a significant dis-agreement be- 
tween these two models and the Yusifov (120041 model; and 
(iii) significant agreement between our observed stellar warp 
and the Yusifov ( I2004> model. In Sect. 12. ll we showed that the 
Yusifov (2004 1 model is based on a rather incomplete catalog, 
thus the closer agreement of our retrieved warp with this model 
was an un-expected result. To solve this confusion however, 
one must go into the details of the models. 

The 3 models describe the warp as a series of concentric 
rings which intersect the Galactic plane along a line of nodes 
(usually measured as the Galactocentric angle from the sun 
to the Galactic center <p, see Fig. and Fig. 1A of Evans et 
al. [R>98 1. The models of Lopez-Corredoira et al. {2002) and 
Drimmel & Spergel ( 12001 \ report <p = -5° and 0° respectively. 
On the other hand, the Yusifov ( 2004 1 model derives a rather 
higher <p = +15°. <p different than zero means that the Sun is 
not located on the line of nodes, instead it is already inside 
one of the two oppositely warped regions. As a consequence, 
when viewed from the Sun, geometrically the warp maximum 
is nearer to us in one hemisphere than it is in the other. Thus, 

4 Note that in this paper we do not use the Drimmel & Spergel 
( 2001 1 model that is based on dust, since this predicts very high warp- 
ing (cf Fig. 2 of Yusifov l2004l . 



if the warp is traced at fixed distances from the Sun and <p + 0, 
then the observed warp is asymmetric. Yusifov (j20Q4) derives a 
positive (f> value and therefore predicts a significantly asymmet- 
ric warp, with a relatively stronger Southern warp maximum 
amplitude. On the contrary, Lopez-Corredoira et al. (f2002Ti de- 
rive a negative <p value which produces an observable asymmet- 
ric warp, with a relatively stronger Northern warp maximum 
amplitude. Lastly, Drimmel & Spergel J2001i derive = 0° 
and therefore their model predicts a symmetric warp; equiva- 
lent Northern/Southern warp maximum amplitudes. 

Another factor that contributes to the appearance of an 
asymmetric warp is due to the chance location of the Northern 
warp just behind the Norma-Cygnus arm (labeled outer in 
Fig. 0- Thus, extinction in the Norma-Cygnus arm and pos- 
sible variations of the extinction curve (due to the penetration 
of gaseous or dusty regions), coupled with the possibility that 
the Sun may not lie on the line of nodes are all factors that con- 
spire in producing a less pronounced apparent Northern warp 
at heliocentric distances of 2.8 and 7.3 kpc. 

To develop further the particular issue of a symmetric and 
asymmetric stellar warp, we now compare the stellar warp as 
traced by the 3 RGB samples at different heliocentric distances. 
Since our results favor rather high and positive <p values we 
will continue the comparison with the Yusifov model as this 
has shown the best fit to our RGB sample at heliocentric dis- 
tance 7.3 kpc. Figure^]shows clearly that the Yusifov model 
provides a satisfactory match with the stellar warp as derived 
by the R Q ~ 2.8 and ~ 7.3 kpc RGB samples. Figure fT2l also 
shows how the Northern warp amplitude is almost half that of 
the Southern warp, for the two distances. 

In principle, the warp derived at R Q ~ 16.6 kpc should not 
be compared with the Yusifov model (Fig. ll3> since that model 
is not applicable at Rqc ^ 14 kpc. Nevertheless, for the de- 
rived warp at R Q ~ 16.6 kpc we find that the Yusifov model 
predicts maximum warp amplitudes that are not very far from 
our derived values; with amplitude differences of 0.25 and 0.15 
kpc above and below the mid-plane, respectively. Moreover, 
the location of the warp maxima are still close to those of the 
derived data-points, indicating that the warp phase angle is still 
in agreement with observations. Figure^] also shows another 
important finding: only at distances of R Q ^ 16.6 kpc can we 
obtain a good symmetry between Northern and Southern warp 
maximum amplitudes; both being around ~ 0.6 kpc above and 
below the mid-plane. The fact that we obtain a symmetric warp 
at R Q ~ 16.6 kpc also indicates that probing the warp at such 
large distances is affected less by the dependence on (p. 

The passage from the stellar warp at R Q ~ 7.3 kpc to that 
at 16.6 kpc shows an out-break in the Northern warp maxi- 
mum by almost 0.45 kpc. This confirms how probing distant 
regions (which are less affected by grand-design structures like 
the Norma-Cygnus arm) reveals the real entity of the Northern 
warp amplitude. On the other hand, the passage to the stellar 
warp at R Q ~ 16.6 kpc shows a limited increase in the Southern 
warp maximum amplitude by ~ 0.2 kpc. This limited increase 
probably reflects the way the Southern gaseous warp shows a 
saturation effect at ~ 0.75 kpc below the plane, that induced 
previous investigations to consider it constant after Rcc = 14.0 
kpc (Burton 1998). 
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Fig. 11. A comparison of our Galactic stellar disk warp height Z - warp (kpc), as obtained from the RGB sample at heliocentric 
distances of 7.3 kpc points with error bars, with models based on stellar tracers lines. Note the excellent agreement between data 
and models for the Southern warp. The disagreement for the Northern warp is explained in the text. 



4.2. Significance of the retrieved stellar warp 

Figures|5] E3 ^2 El an d El a U° w us to conclude the follow- 
ing: 

- The correspondance of the CMa over-density with the 
Southern warp maximum at I ~ 240°: This conclusion is based 
on both our derivation of the stellar warp (at different distances) 
and independently derived models. This establishes that the 
warp maximum is significantly displaced from I =s 270° at CMa 
distances. 

Now that our results indicate that the Galactic warp at R = 
7.3 kpc can bring the mid-plane ~ 3° below the nominal b = 0°, 
star-count comparisons above and below the plane should be 
made as follows: in the Canis Major direction (7 = 240°) star 
counts at CMa core (b = -8°) should be compared with those at 
b =i +2°. This explains the consistent detection of a CMa over- 
density or peculiar signature when comparison fields are taken 
at b - +8° (Martin et al. l2005l and Bellazzini et al. 120041 One 



further piece of evidence which indicates that the CMa detec- 
tion is the recovery of the Galactic warp Southern hemisphere 
maximum is found in the so-called structure A. Along with the 
CMa over-density, structure A is the second over-density found 
in the Martin et al. (2004) analysis. Its location (in the Northern 
hemisphere) makes it fully compatible with being the Northern 
warp maximum, as also found in this paper. 

Yet another overdensity (in Argo) has recently been an- 
nounced (Rocha-Pinto et al. 2005 ). Our results suggest that this 
over-density cannot be directly correlated to the Galactic warp 
maximum (as suggested in Bellazzini et al. l2006t . Indeed, at 
I =i 290° (~Argo center) the Galactic warp amplitude is less 
pronounced than that at / 240°. To consider a possible ori- 
gin of this over-density, we briefly reconsider the analysis of 
Rocha-Pinto et al. ( 12005 \ . In order to increase their accessi- 
bility to low latitude sky areas and reduce the dependency on 
foreground reddening Rocha-Pinto et al. J20051 > use a Galactic 
model to remove the foreground density. The model is idealised 
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Fig. 12. The Galactic stellar disk warp height as derived from 
the RGB samples at R G = 2.8 and 7.3 kpc compared with the 
Yusifov (20049 model. 
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Fig. 13. The Galactic stellar disk warp height as derived from 
the RGB sample at R Q = 16.6 kpc improperly compared with 
the Yusifov 



(I2004> model (not applicable at these distances). 
Note that for Galactocentric distances greater than 14 kpc the 
model imposes a constant Southern stellar warp so as to be 
consistent with the derived parameters of the gaseous warp. 

in that it assumes a cylindrically symmetric density distribution 
about the Galactic center. Figure|9]shows that the Galactic warp 
amplitude at / 290° is not zero, inconsistent with the assump- 
tions of the Rocha-Pinto Galactic model. Although this effect 
remains to be quantified, and although the warp amplitude is 



not extreme at I ~ 290°, the high star density this close to 
the Galactic disk will amplify the residual Galactic component 
when Galactic populations are removed, necessarily affecting 
the quantification of the Argo over-density. 

- The regularity of the Galactic stellar warp out to Rqc ~ 
20 kpc: We have measured the stellar warp at 3 different dis- 
tances and argued that the only visible perturbation away from 
regularity (a "swelling" in the R G = 2.8 and 16.6 kpc rings) can 
be explained as being due to the presence of the Orion-Cygnus 
arm in lines of sight towards the Northern warp and between 
140° < / < 180°. That is, we conclude that a global and reg- 
ular warp signature is traced to Galactocentric distances of at 
least ~ 20 kpc in the anti-center direction. We emphasise this 
finding, since it is a clear conclusion even though there are two 
factors that could have masked such a signature: 1) the Galactic 
flare (see next Section) and 2) the decreasing stellar density in 
the Galaxy's outskirts. 

- The consistency of the stellar warp with the dust and 
gaseous warp: In addition to being a regular large-scale struc- 
ture, the stellar warp is consistent with that of the interstel- 
lar dust and neutral atomic hydrogen. The consistency might 
have been expected given the close physical correlation be- 
tween these 3 components. However, it is when considering 
accretion scenarios that this finding acquires more significance. 
Accretion of companions can be responsible for generating a 
short-lived warp, and can modify differentially the stars-dust- 
gas warp properties [cf Binney (1992)]. In regards to this, it is 
interesting to recall that the numerical simulations by Helmi 
et al. d2003t . cited to explain the almost circular orbits of 
CMa and Mon. Ring, suggest that the accretion event must be 
relatively young (< 1 Gyr) otherwise any coherent structure 
would be dissolved. Developing consistency between dynami- 
cally young events - a few rotation periods - and global regu- 
larity remains to be investigated in detail. 

- The orientation of the Galactic bar and warp: The re- 
sults presented in this paper define a stellar warp with a 
rather high and positive phase-angle. Interestingly, we note 
that a positive warp phase-angle follows the orientation of the 
Galactic bar. The amplitude of the warp phase-angle (+15°) 
in Yusifov's model is very similar to recent determinations of 
the Galactic bar orientation: Freudenreich ( 1998) derive +14°, 
while Babusiaux & Gilmore ( I20051 derive +22° +5.5°. We have 
shown that the pulsar catalog used by Yusifov is biased towards 
the Galactic inner regions (Rqc ^ 8 kpc, see Fig.^. Whereas 
this could have been a shortcoming in his model, it turns out 
that this inner regions weighted catalog may have led to a better 
determination of an important warp parameter; its phase-angle 
at least in the inner Galaxy. Thus, it may be that the Milky 
Way bar and the inner warp-rings share the same orientation 
angle. This point will be further investigated in Momany et al. 
(in prep.). For the time being, we note that our determination 
of a non-zero phase angle frees us from the uncomfortable as- 
sumption of a fortuitous location of the Sun along the line of 
nodes. 
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5. The Galactic flare 

In the present section, and before extending our discussion on 
the Mon. Ring, we derive another necessary Galactic ingredi- 
ent, flaring of the Galactic disk. 

For this analysis we use the Z/, values of the scale height as 
determined in Sec. 2.3 for the 3 M-giants samples at R Q = 2.8, 
7.3 and 16.6 kpc. The flaring of the Galactic disk (viewed as an 
increase of the scale-height with increasing Galactocentric dis- 
tance) is evident in Figure [H] where we trace the radial trend 
of Z;,. Before discussing the details of Figure El we remind 
the reader that: (i) we used a single exponential in fitting the 
vertical density profile, defining a scale-height; (ii) in the fit- 
ting process we excluded regions at low latitudes so as to avoid 
regions of high reddening corrections around the Galactic mid- 
plane; and most importantly (iii) our three 2MASS M-giants 
samples have an upper vertical limit in latitude of b — +/ -20°, 
and this corresponds to different upper Z probing once the dis- 
tances of the 3 samples is accounted for. This also means that 
we have a different weighting of the thin and thick disk stel- 
lar populations for the three M-giant samples. Table [5] reports 
the adopted limits in the vertical density profile fitting, for the 
three samples. Considering that at the solar distance the scale- 
height of the thin disk is ~ 0.3 kpc while that of the thick disk 
is ~ 1 kpc (Wyse & Gilmore 2005}, it is clear that the sample 
at R Q =i 2.8 kpc is weighted more by thin disk populations. On 
the contrary the R Q 7.3 and 16.6 kpc are weighted more by 
thick disk populations. 

Thanks to the use of different lines of sight at fixed he- 
liocentric distances, our data-points cover a wide range of 
Galactocentric distances and thus are reliable in tracing the 
flare all over the Galactic stellar disk (inner and outer). In par- 
ticular we note the following: 

- For Galactocentric distances in the range 4 < Rcc < 15 
kpc, the R e = 7.3 kpc M-giants sample (filled squares), 
shows only small variations and in general displays an 
almost constant scale-height. Interestingly, the last data- 
points clustering around Rqc — 15 kpc show a saturation 
of the scale-height at ~ 0.75 kpc; 

- For Galactocentric distances in the range 8 < Rcc ^ 25 
kpc, the R G = 16.6 kpc M-giants sample (open squares), 
shows a bimodal behavior with an: (i) almost constant 
scale-height within Rcc S 15 kpc; and (ii) increasing 
scale-height for Rcc > 15 kpc. The inner radii with the 
-constant trend show a mean scale-height of the same or- 
der as that derived from the R Q =7.3 kpc M-giants sample. 
By contrast, at large Galactocentric distances with an in- 
creasing scale-height we note that the data-points do not 
converge to a single value, but show a gradual increase 
with distance out to at least Rqc ~ 23 kpc. The general 
trends reflect that seen in Alard (120001 starred symbols) and 
show fair agreement with the Yusifov model. A quantitative 
parametrization of the observed flare and warp will be sub- 
ject of a future paper (Momany et al. in prep.). 

- In contrast to the almost constant scale-height derived from 
the R = 7.3 and 16.6 kpc M-giants samples for the inner 
regions, the sample at R Q - 2.8 kpc (open triangles) shows 
some differences. Between 5 < Rcc ^ 9 kpc, the R = 2.8 



Table 1. The limits within which we fit the vertical scale-height 
for the three M-giants samples. All values are in kpc. 



Distance 


Lower Z limit 


Upper Z limit 


2.8 


0.15 


0.85 


7.3 


0.45 


2.25 


16.6 


0.60 


3.00 



kpc shows a mean scale-height of 0.35 kpc, almost half that 
found for the R Q = 7.3 and 16.6 kpc M-giants samples (0.65 
kpc). 

This difference in scale-height is mainly due to the differ- 
ent weighting of the thin and thick disk populations when 
extracting the M-giants samples at different distances be- 
tween -20° < b < 20°. This translates into different limits 
when performing the scale-height fitting (cf Table 1). For 
the R Q = 2.8 kpc sample the profile fitting was made be- 
tween 0.15 and 0.85 kpc, clearly more weighted towards 
the thin disk regime. Reassuringly, the overall mean scale- 
height of the R Q = 2.8 kpc sample (~ 0.35) is close to typi- 
cal values for that of the thin disk. On the other hand, for the 
R Q = 7.3 and 16.6 kpc samples the profile fitting starts at 
~ 1.5 times the thin disk scale-height. Therefore, the over- 
all mean scale-height of the R = 7.3 and 16.6 kpc samples 
is an intermediate value (0.65 kpc) that reflects the mixing 
of the thin and the thick disk populations and respective 
scale-heights. Lastly, besides the small difference in mean 
scale-height, one notes that the sample at R = 2.8 kpc 
shows a similar trend with Galactocentric distance as does 
the flare model derived by Lopez-Corredoira et al. J2002> . 

In general, our data-points reflect the fact that the scale- 
height increases with increasing Galactocentric radii. For the 
inner Galaxy the scale-height is rather similar to the that in 
the solar neighborhood, a consistency that continues to about 
Rcc ~ 15 kpc. Further out the scale-height grows relatively 
rapidly, with a smooth increase until about ~ 23 kpc. 

The thickening of the stellar disk is the re-distribution of 
stars from lower to higher latitudes. A disk thickening or flar- 
ing is a phenomenon that requires a reduction below single ex- 
ponential extrapolation in the stellar density near the Galactic 
plane for the outer disk. Thus, our detection of (i) a large-scale 
warp signature out to Rcc ~ 24 kpc; and (ii) a disk flaring 
out to Rcc ~ 23 kpc provides new evidence of a rather ex- 
tended Milky Way stellar disk. Therefore, our results add to 
earlier suggestions by Lopez-Corredoira et al. J2002t and Alard 
J2000t that the thin disk cut-off at R C c ~ 15 kpc (Ruphy et al. 
119961 "is not only unnecessary but also inconsistent with the 
data" . 

In Fig. [B] we compare the thickening/flaring of the stel- 
lar disk (as obtained from the RGB samples at R G = 7.3 and 
16.6 kpc) with that found for the neutral hydrogen gas layer 
and molecular clouds ensemble (Wouterloot et al. [1990 1. On 
the one hand, the observations by Wouterloot et al. (1990J are 
a reminder of the presence of molecular clouds with embedded 
star formation out to Rcc ~ 20 kpc. The authors argue that the 
lack of CO emission at Rcc > 20 kpc is not to be attributed to 
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Fig. 14. Variation of the scale-height (derived from the R G = 2.8, 7.3 and 16.6 kpc RGB samples) as a function of Galactocentric 



distances. We also show a comparison of our results with two models and data-points from Alard J2000> . 



the sensitivity of the IRAS survey, rather it shows the absence 
of recent star formation at these distances. On the other hand, 
the Figure proves the compatibility of the stellar and gaseous 
flaring for the outer Milky Way. We remind the reader that the 
filled squares refer to a mixture of thin/thick disk stellar pop- 
ulations, and this explains the gas-stars differences in relative 
scale-heights. The different populations shown have different 
vertical velocity dispersions, hence their different scale heights. 
The only unexpected aspect of this comparison is the very high 
outer gas scale height, apparently exceeding that of the stars 
beyond Rqc ~ 20 kpc. This is not understood. 



6. The warped and flared Galactic disk and the 
connection with stellar over-densities 

Having demonstrated the existence and determined the proper- 
ties of the stellar warp and flare, we now apply this description 
of a deformed stellar disk to investigate its impact on reported 
stellar over-densities at low Galactic latitudes. In the previous 
sections we showed that the warp and flare are fairly well de- 
scribed by the Yusifov model. We therefore will continue us- 
ing this model. We remind the reader that the model describes 
rather well the inner Rcc ~ 14 kpc, but represents only a qual- 
itative description of the outer Milky Way warped and flared 
regions. 
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Fig. 15. A comparison of the thickness of the stellar disk, neu- 
tral hydrogen gas layer and molecular clouds ensemble for the 
outer Milky Way. Filled squares show the flaring for the RGB 
samples at R G = 7.3 and 16.6 kpc, whereas the open squares 
plot the RGB samples at R Q = 2.8 kpc. 



6.1. The Canis Major over-density 

To show tangibly the importance of allowing for deviations 
away from symmetry around b — 0° for studies of the Milky 
Way stellar disk, in Fig. ^] we show a cut in the YZ plane of 
the warped and flared Galaxy 5 using the Yusifov warp model 



See also Fig. 1 A and IB of Evans et al. 119981 
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in the direction [(/, b)= 248°, -7°] and distance [R Q 7.2 kpc] 
of CMa (as most recently derived by Bellazzini et al. l2006> . 

The solid line traces the mean warped stellar disk, whereas 
the grey dashed lines follow the three scale heights of the 
Galactic disk. The Figure shows clearly that the CMa over- 
density falls at less than ~ 500 pc from the Galactic warped 
mid-plane, that is ~ 50% the maximum density of the stellar 
disk at that distance. 

This figure clarifies how important is the Galactic warp to 
any analysis of the CMa over-density. Specifically, comparing 
star counts at Z - 1 kpc with those at Z + 1 kpc one ends 
up comparing a region falling within ~ 50% of the stellar disk 
density with another having only ~ 14%. 



4 

2 

o 

& 



-4 T 



(l,b)=(248;-7') 




-20 



-10 10 

R GC (Kpc) 



20 



Fig. 16. A cut in the YZ plane of the warped and flared Galaxy 
using the Yusifov warp model in the direction and at the dis- 
tance of CMa. The thick line marks the mean warped stellar 
disk whereas the grey dashed lines trace the density at lx, 2x 
and 3x the scale-height of the disk. 



6.2. The Monocerous Ring over-densities 

The knowledge of the outer structure of the Galactic disk de- 
rived in this analysis does not provide an explanation for the 
Monocerous Ring. However, we show that the Galactic flare 
can account for some aspects of some positive detections of 
the Mon. Ring in various surveys and therefore, needs to be 
taken into account in investigations of the relative importances 
of accretion or Galactic structure for the origin for the Mon. 
Ring. 

In Fig. ^] we show the location of positive Mon. Ring de- 
tections with respect to the warped and flared Milky Way disk. 
For the Sloan fields (Newberg et al. 120021 1. we plot detections 
at heliocentric distances of 7, 8, 9 and 10 kpc, whereas for the 
INT fields (Conn et al. l2005l and Ibata et al. l2003l we plot the 
detections at the specific distances given in Table 3 of Conn et 
al. 

Figure ^] shows that all reported Mon. Ring detections 
are at locations where the structure of the outer warp and 
flare contribute to the analysis, and must be carefully consid- 
ered, and at locations where the dynamical interaction between 
Monocerous and the outer disk will be very important. The 
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Fig. 17. The Same as in Fig ^] but for fields where positive 
detections of the Mon. Ring have been reported. For the Sloan 
detections (upper three panels) we plot detections at heliocen- 
tric distances of 7, 8, 9 and 10 kpc. Similarly, the lower three 
panels show INT detections (Conn et al. 2005 and Ibata et al. 
120031 plotted at the distances given in Table 3 of Conn et al. 



Mon. Ring has attributed properties (a radius of 15 - 20 kpc, 
a radial thickness of 2 kpc and a vertical scale-height of 0.75 
kpc, Ibata et al. 12003b which place it well within the outer disk 
defined earlier in this paper. 

Of all the Mon. Ring fields, the one at {I, b) = (150°, +15°) 
is perhaps the most important. In their INT survey, Conn et al. 
(2005) argued that the detection of Mon. Ring structures on 
both sides of the Galactic plane, at similar Galactic longitudes, 
is the strongest evidence that the Mon. Ring is not an artefact 
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of Galactic warp effects on star counts. We agree with Conn 
et al. (©05 1 that inadequate substraction of the warp is not a 
valid explanation for Mon. Ring detections above and below 
the mid-plane. We consider here if disk flaring provides a valid 
explanation for the Mon. Ring detections in both hemispheres 
at (l,b) = (150°, ±15°). 

A hint that Mon. Ring detections are related to disk flaring 
comes from the field at (l,b) = (123°, -19°), first studied by 
Ibata et al. J2003I I. who reported a positive Mon. Ring detection 
at Rqc ~ 14 kpc. This field was re-examined by Conn et al. 
J2005> who reported a possible new detection at Rcc ~ 21 kpc. 
Figure^]shows that the earlier and this new tentative detection 
of Mon. Ring features follow the expected Z height of the flared 
MWdisk. 

We are clearly far from performing a quantitative analysis 
of the impact of a flared disk on positive detections of Mon. 
Ring. On the one hand, we should do our best to decrease the 
rather large uncertainties in the distances of Mon. Ring detec- 
tions. On the other hand, it would be interesting to have new 
Galactic models in which, for example, the radial truncation 
of the stellar disk at Rqc ~ 14 kpc is dropped, the orientation 
of the warp phase-angle is updated, and new flare modeling 
is included. Note for example that the extensive comparison 
between observed and Besancon simulated diagrams in Mon. 
Ring fields (Conn et al. 120051 implicitly assume a stellar disk 
truncation at Rqc ~ 14 kpc. Since most of the Mon. Ring de- 
tections are at Galactocentric > 14 kpc it follows that the com- 
parison with the Galactic model is inconsistent. 

7. Galactic structure to R GC = 20 kpc in the outer 
Galaxy: what do we expect? 

In the previous Sections we have shown how complex the outer 
(beyond the Solar circle) stellar disk of our Galaxy is, and how 
many intrinsic structures can be identified in it. We have also 
critically discussed the difficulties in determining the stellar 
disk main parameters (warp, flare, scale height) by studying 
its main stellar components. Despite these unavoidable diffi- 
culties, we have shown that the parameters inferred from the 
stellar distribution compare well with what has been obtained 
from a number of investigations on the distribution of the gas, 
dust, and neutral Hydrogen. In order to complete our discus- 
sion, in this section we add additional information on what we 
know about the outer Galactic disk, and some information of 
the properties of the stellar disk in other galaxies. 

The outer structure of galactic disks, and the Galactic disk, 
has been extensively studied for very many years, prior to 
the current re-discovery of its intrinsic interest. A consider- 
able amount of direct information on the outer Galaxy is avail- 
able, as is further information by analogy from studies of other 
galaxies. 

A superb review by Burton ( 119881 1 summarizes and ex- 
plains the subject and data, while presenting the observed ba- 
sic features and results which remain valid in the light of new 
data. The most extensively studied major component is the H i. 
This H i extends approximately smoothly and continuously to 
R ~ 25 kpc (Burton [THHH Figs 7.19, 7.20). The distribution 
is significantly warped and flared (Figs 7.18 to 7.23), there 



are significant non-circular motions, and the gas is distributed 
asymmetrically in azimuth (Fig 7.21), implying a very sig- 
nificantly lop-sided gas distribution. The results from Burton 
J1988I I were updated and extended by Hartmann & Burton 
i 19971 1. and most recently in a very extensive new study of the 
whole sky Hi distribution, using significantly improved cali- 
brations and stray-light corrections, by Kalberla et al. J2005t . 
Inspection of the Hartmann/Burton maps for velocities corre- 
sponding to the far outer Galaxy (Vlsr ^ +50 km/s) shows 
extremely clearly that the outer gas layer exists, is smooth in 
velocity space (i.e. distance), is centered significantly below the 
Galactic plane, is lop-sided and warped, and is thick. By con- 
trast, until the present resurgence of interest, there have been 
rather few large angular-scale studies of the stellar distribution. 

A very important point to note in anti-center sight-lines is 
the considerable 'foreground' complexity associated with the 
spiral arm structure, which extends well beyond the solar cir- 
cle: in particular, there is the Perseus spiral arm just outside the 
solar orbit, clearly visible in CO / - v maps extending across 
the entire anti-center region from / « 60° through 180° to 
/ w 260°. In the anti-center foreground, up to 2 kpc distant, 
especially in the directions towards I = 140° - 220°, there 
are large angular-scale star-forming regions, with their asso- 
ciated dust, including Monoceros, Orion, the Taurus-Perseus- 
Auriga complex, and UMa and Cam OB-star regions. Some of 
these large-angle structures are physically large in space (e.g. 
Gould's belt) and are distributed highly asymmetrically with 
respect to the main Galactic Plane. 

Thus the simplest and most robust conclusion from a direct 
examination of the sky is that one does not expect any simple 
symmetric distribution of sources on the sky to be an adequate 
description of the outer Galaxy. One must select tracers by dis- 
tance, and appreciate complex foregrounds, and expect the out- 
ward continuation of normal spiral structure. 

7.1. Hi and CO surveys: a star-forming disk to 
R GC = 20 kpc 

Diplas & Savage (I199U . in an early extensive study of outer- 
Galaxy gas which remains valid and topical, conclude: for 
12 < Rcc < 18 kpc the z-distribution of the Hi is complicated, 
and reveals the existence of a confined component and an ex- 
tended component. The confined component exists only within 
a Galactocentric radius Rcc ~ 18-20 kpc, while the extended 
component can be followed in some directions to Rcc ~ 25-30 
kpc. In the directions of maximum warp (90° < I < 110°) the 
average distance of the gas away from the the Galactic plane is 
found to reach ~ 4 kpc at Rcc ~ 24 kpc, and the flaring of the 
layer as measured by z rms increases to ~ 3 kpc. In these direc- 
tions the warping outer layer appears to be reliably traced over 
a region at least 10 kpc thick in z. When viewing from the po- 
sition of the Sun, this implies that the warped and flaring outer 
Galaxy extends to Galactic latitudes as large as 25°. 
This raises the interesting question: is the "confined compo- 
nent" what has been re-discovered recently, the Mon. Ring ? 

Nakanishi & Sofue d2004t completed a global analysis of 
Galactic H i: they show extensions of H i in the outer Galaxy 



18 



Y. Momany et al.: The Galactic stellar warp and flare 



well beyond 20 kpc, with evidence of very large scale rotational 
asymmetry (lop-sided), so much so that the Galaxy extends far- 
ther in the first and fourth quadrants. An important point rele- 
vant to current discussions is that the general outer distribution 
of the H i is connected but is not simply smooth. Since the dis- 
tances are derived from a kinematic model (flat rotation curve), 
there is a direct degeneracy between kinematic asymmetries - 
i.e. elliptical orbits, and spatial inhomogeneities. Outer spiral 
structure is evident, as is vertical warping. 

We emphasize that from all these and the many other Hi 
and CO studies there is no evidence for a significant perturba- 
tion from a major satellite accretion event currently underway: 
quite the opposite in fact, the global approximate smoothness 
of the H i distribution is very difficult to reconcile with such a 
model. 

There are many molecular studies of the outer Galaxy indi- 
cating that the gas undergoes "normal" outer galaxy rates and 
distributions of star formation, again supportive of a largely 
steady-state system. An example is the Brand & Wouterloot 
(1995) study of outer IRAS sources, which are star-forming 
molecular clouds out to Rqc — 20 kpc. They collect data from 
several earlier studies. Snell et al. (I2002> provide a valuable 
overview of both the Galactic data and studies of other galax- 
ies, concluding that outer Galaxy molecular cloud properties 
and star formation is a natural extension of the inner Galaxy 
properties. 

We note a study by Freudenreich (J1996J of the DIRBE ob- 
servations of outer galaxy dust, which reports a detection only 
to Rqc - 12 kpc. It is unclear why this study fails to detect 
the farther Galactic molecular material, but this paper is the 
anomaly, all the other studies agree rather well. 

A major CO study by Dame, Hartman & Thaddeus ( 120011 
and recent extension to higher latitudes by Dame & Thaddeus 
( 2004) shows extended CO emission at latitudes up to 50° in 
the anti-center region, with this gas probably being associated 
with star formation in Taurus-Auriga. Their Fig. 3 shows very 
well the outer molecular gas at velocities up to + 100 km/s, cor- 
responding to very large Galactocentric distances. 

The most recent high-sensitivity southern Hi survey, 
Dickey et al. (2004i shows that the Hi warp continues to more- 
negative z-distance as distance increases to a minimum be- 
tween Rcc = 13-16 kpc, with the location a function of 
longitude through this distance range, and then returns to near 
Z = by Rqc = 20 kpc. Beyond 20 kpc and out to 30 kpc 
the Hi distribution 'is dominated' by a density enhancement 
which extends coherently over very wide angles, and may be 
a coherent extension of the Perseus spiral arm for longitudes 
I « 220° - 325°. This very large structure, if indeed it is coher- 
ent, may be similar to the spiral-like structures in H i extending 
well beyond the visible optical disk seen in several other galax- 
ies, with associated star formation and apparently normal struc- 
tures, by Ferguson et al. dl998b and the several other examples 
recently observed by UV analysis. 

As an example of studies of external spirals, we note 
the complex, asymmetrical, inflowing Hi seen in NGC2403. 
Fraternelli et al. d2002> show this is a normal, non-interacting 
spiral galaxy, suggesting that such H i features may be common 



among spiral galaxies and that, perhaps, they have not been de- 
tected yet because of the low sensitivity of previous surveys. 

7.2. Expected kinematics 

Bailin & Steinmetz ( 120031 1 note that warps are very common, so 
much so that perhaps 50% of galaxies have them (Reshetnikov 
& Combes (119981 . But, warps are not long-lived (e.g. Binney 
I1992> . so they must be generated in some way. 

Hence, if a warp is kinematically driven, there is no reason 
to expect simple geometry or kinematics. The actual kinematic 
pattern will depend on what is the cause, which is not clear. 
Possible causes include the cumulative effects of many satel- 
lites (not necessarily on decaying orbits), though in the case of 
the Galaxy it may just be the LMC which induces a large-scale 
reaction from the DM halo (Weinberg [T9981 Tsuchiva 120021 
but see Garcfa-Ruiz, Kuijken & Dubinski 2002) or maybe by 
Sagittarius (Bailin & Steinmetz l2003> . Bailin & Steinmetz em- 
phasize that a warp generated by interaction with a satellite is 
a dynamical warp, so that one does not expect to see stationary 
or uniformly precessing kinematics. A transient induced warp 
moves outward at the group velocity of bending waves in a stel- 
lar disk, which scales as n G X E/f2, for surface density £ and 
rotational velocity Q. (Hofner & Sparke [T994t or ~ 30kpc/Gyr 
near the Sun. 

Newberg & Yanny ( 2005) suggest significant triaxiality in 
the large scale structure of the outer galaxy star counts. This is 
for distances of 3 - 16 kpc exterior to the Sun. This follows ear- 
lier suggestions by Roberta Humphreys. Note that they see sim- 
ilar triaxiality and coherent radial motions in the local old disk, 
a pattern perhaps consistent with the combined perturbations 
on circular orbits due to spiral structure (cf de Simone, Wu & 
Tremaine 2004it and effects due to a bar-shaped perturbation in 
the disk/bulge potential (cf Babusiaux & Gilmore 2005 1. In any 
event, simple circular orbits are neither natural nor expected. 

7.3. Expectations from other galaxies 

Star formation in very outer disks, beyond two R25, has 
been detected from H-alpha imaging in several galaxies (e.g. 
NGC6 28, NG C1058, NGC6946 Ferguson et al. 119951 . Thilker 
et al. ( 120051 GALEX imaging the outer disk of M83, shows 
current star formation at up to three RHolmberg m the warped 
outer disk of M83. The extreme patchiness of the very outer 
disk is evident in Fig. 1 of Thilker et al., yet there is no sug- 
gestion that what is observed is anything other than star for- 
mation in a disk, albeit at very large Galactocentric distances. 
The GALEX images show star formation, and the existence of 
both young and intermediate-age open clusters, well beyond 
an apparent disk 'edge'. The authors note the location of the 
young star-forming regions "preferentially on local maxima or 
filaments in the structure of the warped Hi disk", with these 
complexes superimposed on an extended field stellar popula- 
tion. 

Gentile et al. d2"0"0"3l looked at ES0123-G23, an apparently 
thin edge-on galaxy with a warp in the optical, and saw Hi 
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warped gas up to 15 kpc from the optical major axis, implying 
a warp through 30°. 

An interesting specific study, which reminds us that the dy- 
namics of outer galaxies may not always require current ac- 
tive mergers, is that by Fraternali et al. ( 120021 of NGC2403. 
They studied the spiral galaxy NGC2403 with the VLA, and 
discussed the properties of the extended, differentially rotating 
Hi layer with its Hi holes, spiral structure, and outer warp. 
In addition, their new data reveal the presence of a faint, ex- 
tended, and kinematically anomalous component. This shows 
up in the H i line profiles as extended wings of emission toward 
the systemic velocity. In the central regions these wings are 
very broad (up to 150 km/s) and indicate large deviations from 
circular motion. They separated the anomalous gas component 
from the cold disk and obtained a separate velocity field and 
a separate rotation curve for each. The mass of the anomalous 
component is 1/10 of the total Hi mass. The rotation velocity 
of the anomalous gas is 25-50 km/s lower than that of the disk. 
Its velocity field has non-orthogonal major and minor axes that 
they interpret as due to an overall inflow motion of 10-20 
km/s toward the center of the galaxy. The picture emerging 
from these observations is that of a cold H i disk surrounded by 
a thick and clumpy H i layer characterized by slower rotation 
and inflow motion toward the center. The origin of this anoma- 
lous gas layer is unclear. It is likely, however, that it is related 
to the high rate of star formation in the disk of NGC2403 and 
that its kinematics is the result of a galactic fountain-type of 
mechanism. 

In a related study, Fraternali et al. d2004> show that thin 
'cold' Hi disks are often surrounded by thick layers (halos) of 
neutral gas with anomalous kinematics. They present results for 
three galaxies viewed at different inclination angles: NGC891 
(edge-on), NGC2403 (i = 60), and NGC6946 (almost face-on). 
These studies show the presence of halo gas up to distances 
of 10 - 15 kpc from the plane. Such gas has a mean rotation 
25 - 50 km/s lower than that of the gas in the plane, and some 
complexes are detected at very high velocities, up to 200 - 300 
km/s. The nature and origin of this halo gas are poorly under- 
stood. It can either be the result of a galactic fountain or of 
accretion from the intergalactic medium. It is probably analo- 
gous to some of the High Velocity Clouds (HVCs) of the Milky 
Way. 

The overall picture from these many studies is that the outer 
parts of galaxies are complex, asymmetric and time-dependent. 
While not understood, ongoing mergers of significant compan- 
ions seem not to be the primary explanation or cause at work. 

8. Discussion and Conclusions 

We have used the 2MASS all-sky catalog to probe two impor- 
tant outer Milky Way parameters; the Galactic warp and flare. 
We show that the Galactic warp can fully explain the detec- 
tion of the Canis Major over-density. We also show that the 
Galactic flare of the outer disk will be apparent at those places 
where there are detection of Monocerous Ring over-densities. 
Thus, the results presented in this paper provide a reminder 
that disentangling the accretion history of the outer disk from 
its complex and time-dependent dynamics remains a challenge. 



Our results add to the well-known evidence of the complexity 
of outer Milky Way structure. 

The Galactic warp: using 2MASS red clump and red gi- 
ant stars we first demonstrated the dependency of the warp 
parametrization (amplitude and phase-angle) on the contamina- 
tion status of the adopted stellar tracer. Red giants are proved 
excellent tools in probing the outer most parts of the Milky 
Way. The detection of a global warp signature by RGB stars at 
mean heliocentric distances of R Q = 2.8, 7.3 and 16.6 kpc, has 
helped us understand the following: 

- The warp phase-angle is oriented in such a way that at fixed 
distances from the Sun we observe an asymmetric warp 
(relatively stronger maximum amplitude in one hemisphere 
with respect to the other hemisphere). For the R Q = 2.8, 
7.3 kpc RGB samples, the Southern hemisphere warp max- 
imum shows twice the amplitude of that in the Northern 
hemisphere. A symmetry of the warp maximum is obtained 
only for the R e = 16.6 kpc, at which the dependency on the 
phase-angle is weakened; 

- The warp phase-angle (</> =* +15°) is close to the orientation 
of the Galactic bar; 

- The detection of the stellar warp in the RGB sample at a 
mean distance R Q = 2.8 kpc proves that the warp starts 
within the solar circle; 

- The detection of the stellar warp in the RGB sample at a 
mean distance R G = 16.6 kpc proves the presence of an ex- 
tended Milky Way stellar population out to Rcc ~ 24 kpc. 
This is in clear contradiction with the assumption of a ra- 
dial truncation of the stellar disk at Rcc ~ 14 kpc. The de- 
tection of an extended stellar disk is confirmed by parallel, 
and independent studies reporting: (i) recent star formation, 
traced by molecular clouds, out to Rcc ~ 20 kpc; and (ii) an 
extended H i spiral arm in the Third and Fourth quadrants; 

- For the 3 RGB samples, the Southern hemisphere warp 
maximum is found around / = 240°, we therefore establish 
the identity of the CMa center with the warp maximum. 

- The derived stellar warp shows excellent agreement (in 
both the phase-angle and amplitude of the maximum) with 
that of the gaseous and interstellar dust distribution. This 
resolves conflicting evidence of amplitude differences be- 
tween the 3 Galactic components. The origin of the differ- 
ent warp signatures is due to foreground contamination and 
distance resolution of the adopted tracer; 

These results do not answer the fundamental question of 
what is the driving mechanism of the warp ? There is no gen- 
eral agreement on which process, or processes, generate warps 
(see Castro-Rodriguez et al. 12002b . Satellite accretion remains 
one among many scenarios put forward to explain warps. Other 
possibilities include: cosmic infall, bending instabilities, in- 
tergalactic magnetic field, accretion of intergalactic medium 
(see also Olano d2004> for a possible warp-HVCs-Magellanic 
Clouds connection). Any proposed scenario in a special case 
should also account for the fact that warps are a ubiquitous 
property of most, perhaps all, disk galaxies (Sanchez-Saavedra 
et al. l20l)3l . 

The Canis Major Overdensity: We have shown here that 
CMa is the southern hemisphere maximum stellar warp. We 
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also note that all CMa properties (measured so far) point to the 
over-density reflecting a normal warped disk population rather 
than that of a cannibalized dwarf galaxy. Among these we note: 

- the absence of any peculiar radial velocity signature (see 
Appendix. A); 

- a possible shortage of an old (>9 Gyr) and metal poor stel- 
lar population (see Appendix. B); 

- the absence of any peculiar proper motion signature (see 
Appendix. C) and a rather problematic orbital motion 6 ; and 

- a rather high line-of-sight extent (see Appendix. D). 

The Galactic flare: having characterized the warp, we also 
derived the Galactic flare for the 3 RGB samples. Unable to 
kinematically discern thick from thin disk populations, the re- 
trieved variation with Rcc of the disk thickness shows a rather 
constant scale-height (~ 0.65 kpc) within Rcc ~ 15 kpc. 
Further out, the mean scale-height increases gradually reach- 
ing ~ 1 .5 kpc at Rqc ~ 23 kpc. Whereas this trend refers to a 
mixture of thin-thick disk populations, for the RGB sample at 
R Q ~ 2.8 kpc (dominated by thin disk populations) we derive a 
mean scale-height of ~ 0.35 kpc that is consistent with typical 
values for the thin disk. 

Overall, we trace the stellar disk flaring out to Rcc ~ 23 
kpc. Thus, we add further and final evidence on the inconsis- 
tency of the disk radial truncation hypothesis, at Rcc ~ 14 kpc. 
Indeed, the detection of stellar disk warping and flaring out to 
Rcc ~ 23 kpc is the key answer to the observed drop in stellar 
density near the Galactic plane. 
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Appendix A: The absence of peculiar CMa radial 
velocity signature 

In "Why the Canis Major over-density is not due to the Warp", 
Martin et al. (2004) presented radial velocities of M-giants at 
the center of CMa. The overall distribution of radial veloci- 
ties was strongly characterized by a double-peak morphology. 
Martin et al. used this distribution to argue that the peak at 
V m d = 109 + 4 km/s is to be attributed to an extra population, 
CMa, that is distinguishable from typical disk stars peaking at 
V rad = 61 +4 km/s. 

A recent analysis by the same group (Martin et al. 12005b of 
more than 1500 RGB and RC stars, showed that the bimodal ra- 
dial velocity distribution was artificially produced by template 
issues. The new V ra d distribution is now broad; there is no pe- 
culiar signature of CMa stars and these are indistinguishable 
from surrounding disk stars. Nevertheless, Martin et al. pro- 
ceeded and have identified a "peculiar distance- radial veloc- 
ity relation that is unlike that expected from thin or thick disk 
stars". That is, while CMa stars show a mean radial velocity of 
72 + 7 km/s at R =5.5 kpc, the mean V ra d increases to about 
114 + 2 km/s at R Q =8.5 kpc. 

To date, this "peculiar" R Q - V ra d relation is considered the 
strongest kinematic evidence standing in favor of an accretion 
scenario. However, as we shall demonstrate below, this relation 
reflects nothing more than the Galactic differential rotation. 

To shed more light on this "peculiar" R Q - V m d relation, 
we combine the results obtained by Martin et al. (their Table 
2) with a kinematic simulation from the Besangon Galactic 
model. In the left panel of Figure IA.H we show the R Q - 
V m d distribution for all simulated stars in the simulated field 
centered at (/, 6)=(240 o ,-8°). We remind the reader that the 
Besangon model imposes a stellar disk truncation already at 
R Q = 7.7 kpc (R cc = 14.0 kpc) 7 , more or less at the CMa 
distance. Thus, at about R e = 7.7 kpc and beyond, one can- 
not quantitatively compare the observed kinematics with the 

7 As demonstrated in Section|5| the results presented in this paper 
confirm earlier suggestions that a disk truncation hypothesis at Roc — 
14.0 is neither needed nor justified. 
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Fig. A.l. Explaining the "peculiar" R e - V ra d relation found for CMa stars as due to Galactic differential rotation. The left panel 
shows a Besancon simulated R e - V m i distribution for stars at the CMa center. Grey dots highlight the modeled Thin Disk 
truncation at R = 7.7 kpc (marked by a horizontal line). Black dots plot the Thick Disk and Halo stars. Open triangles show the 
R Q - V ra d distribuition for simulated stars extracted from the synthetic color-magnitude diagram (upper right panel), and falling 
within a box delimiting the RC region between 5 < 7? < 8 kpc. The open squares are the observed-RC sample from Martin et al. 
(their Figure 6). Right middle panel shows the distribution of radial velocities of the Martin et al. observed-RC sample (dashed 
histogram), and the scaled distribution of simulated-RC stars (continuous histogram). The right lower panel shows the mean V ra d 
of the dwarfs-decontaminated and RC-simulated sample (open triangles) estimated in three distance intervals. The open squares 
plot the observed mean V ra d values for CMa RC stars (see Table 2 of Martin et al.) in the same distance intervals. The error bars 
are from the same table. The Dashed line is a linear fit to the Martin et al. mean values; which is fully compatible with recent 
estimates of the A Oort constant being A = 17.23 + 0.35. 



incomplete Besancon simulation. To highlight the abrupt trun- 
cation of the simulated Galactic disk, Thin Disk stars within 
R Q — 7.7 kpc are plotted as grey dots, while those further out 
(Thick Disk and Halo) are plotted as black dots. Before going 
into more detail, and prior to any selection or analysis, one al- 
ready sees that there is an intuitable R Q - V ra d trend in the sim- 
ulated data, similar to that found in the Martin et al. analysis 
for CMa stars. 

The RC sample of Martin et al. (their Fig. 6) is presented in 
the left panel (open squares). In order to compare this observed- 
RC sample with an appropriate counterpart from the Besancon 
simulation (a simulated-RC sample), we extract RC stars from 
the one square degree, dereddened K a , (J - K) a diagram. The 
simulated color-magnitude diagram is displayed in the upper 



right panel. In grey symbols we highlight all simulated stel- 
lar populations between 5 < R G < 8 kpc. The selection box, 
where we extract simulated-RC stars between 5 < R < 8 
kpc is shown as grey rectangle area. We emphasize that this 
simulated-RC sample, as that studied by Martin et al., is sub- 
ject to foreground contamination. 

The advantage of using simulated data is that one can easily 
isolate the foreground dwarf contamination from the simulated- 
RC sample, by means of the simulated stars distance entries. 
All stars included in the simulated-RC sample are plotted in 
the left panel (open triangles), showing clearly the contamina- 
tion by local dwarfs within ~ 1 kpc. Between 5.5 and 8.0 kpc 
simulated-RC stars overlap with the Martin et al. observed sam- 
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pie (a colored version of the figure is more appropriate for the 
disentangling). 

The middle right panel shows the radial velocities distri- 
bution of the Martin et al. observed-RC sample (dashed his- 
togram). The distribution of simulated-RC stars as extracted 
from the Besancon simulation is plotted as a continuous his- 
togram. We make use of the Martin et al. J2005l > estimate of the 
V m d intrinsic dispersion of disk stars at the CMa distance and 
add this value (via a Gaussian distribution) to the simulated-RC 
stars. The model histogram has been scaled to fit the observed 
histogram: other appropriate scaling is not possible since the 
Martin et al. ( 2005 1 selection function of RC stars is unknown. 
The similarity of the two distributions is remarkable, we note 
however that the histogram of simulated-RC stars shows a 
slightly higher distribution (with respect to the observed sam- 
ple) for stars with V ra d < 80 km/s, and a lower distribution for 
Vrad > 80 km/s. The later feature can be explained as being due 
to the truncation of disk population in the Besancon model at 
R Q = 7.7 kpc: simulating more distant populations will result in 
more stars only having V T ad higher than 80 km/s. On the other 
hand, the higher distribution of simulated-RC stars at V ra d < 80 
km/s can be due to a stronger foreground dwarf contamination 
in the input model (note in fact the discrete and distinguishable 
distribution of local dwarfs within 1 kpc). The two distributions 
remain, however, very similar and this is further demonstrated 
in the lower right panel where we (i) exclude nearby dwarfs 
from the simulated-RC sample, (ii) use this cleaned simulated- 
RC sample to estimate the mean V ra d in 3 distances intervals 
[5 < R Q < 6, 6 < R B < 7 and 7 < R & < 8 kpc]; and 
(iii) compare the cleaned and simulated V ra d values with those 
derived by Martin et al J2005> . in the 3 distance intervals. 

We emphasize the fact that the mean CMa V ra d values de- 
rived in the 3 distance intervals (Martin et al 2005 open squares 
in the lower right panel) are based on what they call "CMa 
RC sample", which they disentangle from the "contaminating 
population". Keeping in mind all the uncertainties associated 
with the Besancon model, one clearly sees that the cleaned 
simulated-RC sample (open triangles in the lower right panel) 
shows, within Icr, the same R Q - V ra d trend found by Martin et 
al. for a cleaned CMa RC sample. 

One last piece of evidence which demonstrates that the 
"peculiar CMa observed R - V m d trend" is nothing but the 
imprint of Galactic differential rotation of normal disk stars 
is shown in the lower panel of Fig. IA.1I Indeed, a linear fit 
of the three points reported by Martin et al. ( I2005l > yields 
R Q = 0.067 x V ra d- Using the Oort equation: 

V rad = A x R B x sin(2 x 1°) (A.l) 

V rad -Ax sin(2 x 240°) x R = (1/0.067) x R Q (A.2) 

one easily gets A = 17.23 + 0.35, fully compatible with 
recent estimates of the A Oort constant (see Dehnen 1998 1. 
Obviously, one has to consider the CMa large heliocentric dis- 
tance and location with respect to the Galactic disk, however, 
the linear relation is really surprising, and does not leave room 
for any other easy interpretation. 

Given that (i) the bimodal V ra d distribution (Martin et al. 
I2004t was due to an artifact in their reduction procedure, and 



(ii) that the claimed peculiar R Q - V ra d relation (Martin et 
al. 120051 is well reproduced by a self-consistent kinematical 
model, and most importantly, seems to be nothing more than 
the imprint of differential Galactic rotation. We therefore con- 
clude that the stellar populations of CMa do not show any pe- 
culiar or distinguishable V ra d signature with respect to normal 
disk stars. 

Appendix B: The possible shortage of an old (>9 
Gyr) and metal poor CMa stellar population 

All studied dwarf galaxies are known to possess an ancient stel- 
lar population that is traced by the detection of either > 9 Gyr 
RR Lyrae variables and/or red/blue horizontal branch stars of 
similar age. Dwarf galaxies in the Local Group show no ex- 
ception to this rule [Mateo[l998| Tolstov l2000l Grebel l200Tl 8 . 
To this adds the most recent entry, Ursa Major (Willman et al. 
2005 ), which shows a clear clump of blue horizontal branch 
stars. Those in favor of an extra-Galactic accretion scenario 
of CMa often highlight the similarities (e.g. mean metallicity, 
mean mass, My etc) between the CMa stellar populations and 
the Sagittarius dwarf spheroidal. As a consequence, it is partic- 
ularly important to investigate the presence of an old (> 9 Gyr) 
and metal-poor population of CMa. 

Needless to say, an old and metal-poor population is nec- 
essary to metal enrich successive generations of stars, other- 
wise one cannot explain how a metal-rich galaxy ([Fe/H]> -1) 
could have built metal-poor clusters. The Sagittarius dwarf ful- 
fills this requirement; i.e. besides being on average a metal-rich 
galaxy, it possesses a metal-poor component that is traced by 
RR Lyraes and blue horizontal branch stars (see Monaco et al. 
2003 and references therein), and possibly extreme horizontal 
branch stars (Momany et al. l2004at . This explains the progeny 
of metal-poor clusters like Terzan 8 and M54. Studying the 
Sagittarius stellar populations outside the tidal radius of M54, 
Monaco et al. J2003t proved that the group of blue horizon- 
tal branch stars identified in S341+57-22.5 by Newberg et al. 
(20021, stre fully compatible with the blue horizontal branch 
of Sagittarius. This shows that a stream, narrow by definition, 
originating from a disrupting dwarf is likely to leave a footprint 
of all its stellar populations (both red and blue). 

Despite the fact that contamination can seriously hamper 
the identification of old blue horizontal branch stars near the 
Galactic plane, we note that the CMa over-density has not 
been associated with a clear identification, or even a hint, of 
old blue horizontal branch stars. Figure |B~TI displays a wide- 
field color-magnitude diagram of the CMa center, based on 
archival ESO/2.2m telescope data. The data has been reduced 
and calibrated following the standard techniques presented in 
(Momany et al. 1200 It . Upon this color-magnitude diagram we 
over-plot the mean location of red clump and blue horizontal 
branch stars in the Sagittarius dwarf (as derived in Monaco et 

8 Izotov & Thuan <2004> proposed a scenario in which young (~ 
500 Myr) galaxies like IZw 18 are still being born in today's universe. 
Nevertheless, in Momany et al. 1 2005 I we showed that I Zw 1 8 can be 
much older and that, at the moment, there is no photometric evidence 
of the existence of young galaxies. 
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Fig. B.l. A 1° dereddened color-magnitude diagrams of the 
core of CMa. Over-plotted is a 6.3 Gyr Z = 0.004 ([Fe/H]^ 
-0.76) isochrone from The Girardi et al. (2003 library. The 
assumed distance modulus is (m - M) a = 14.5. Each star 
was dereddened according to its position in the Schlegel et al. 
d!998> reddening maps. Also plotted are the red clump and 
blue horizontal branch location of the Sagittarius, as derived 
by Monaco et al. ( I20031 . 



al. l2003> . shifted to match the CMa red clump population. No 
hint of an over-density of stars at the expected blue horizontal 
branch location can be seen. Instead, the distribution of stars 
between 14.0 < V a < 16.5 and (V - I) a < 0.3 smoothly in- 
creases towards fainter magnitudes. We remind the reader that 
at V - 16.0, and in between 0.0 < (V - 7) < 0.5 there lies 
a population of young < 100 Myr main sequences stars which 
belong to the Norma-Cygnus spiral arm, as demonstrated by 
Carraro et al. 120051 by means of UBV two-color diagrams. 
This population is intermediate between the disk MS stars 
[oblique sequence extending from (V - I) - 0.2 at V 12.0 
to (V - T) — 0.6 at V 18.0] and the blue plume population 
at bluer colors. An eventual over-density due to the presence of 
blue horizontal branch star (7, 000 < T eS < 11, 000K) should 
take place in between the disk main sequence and blue plume, 
but this region lacks any significant over-density. A similar in- 
dication is found in the optical color-magnitude diagrams of 
Martmez-Delgado ( 120051 and Bellazzini d2006l . 

Recently, Mateu et al. ( 120046 reported on a preliminary RR 
Lyrae search in an area of 8.3 square degrees at the center of 
the CMa. Five RR Lyraes were found at heliocentric distances 
of less than 6.5 kpc. Assuming the most recent CMa distance 
range (6.2-8.2 kpc by Bellazzini et al. I2006l > the 5 variables 
with a mean distance of 5.6 kpc are outside the main body 
of the CMa over-density. Yet, as noted by the Mateu et al., 
the detection of 5 variables is higher than what would be ex- 
pected in the same volume of the Galactic Halo (~1 variable). 
Unfortunately, a direct comparison with the expected number 



of Disk RR Lyrae is difficult since the distribution (and density) 
of RR Lyrae stars in the thick disk is not known. Moreover, one 
should add the current uncertainty in (i) disentangling Disk and 
Halo RR Lyrae at low Galactic latitudes, and (ii) the disk den- 
sity profile as a function of Rcc- As demonstrated thoroughly 
in this paper, the disk radial extension does not show an abrupt 
truncation at Rcc ~ 14 kpc, and this might increase the number 
of expected variable stars. 

We therefore conclude that there is no evidence which 
demonstrates the presence of an old, blue horizontal branch 
population CMa. At the same time, RR Lyrae surveys leave this 
issue rather 2open, highlighting intrinsic difficulties (expected 
Galactic contamination and relative vicinity of CMa) of similar 
studies. More data is needed to draw a firmer conclusion. 

Appendix C: The negative vertical velocity of the 
Canis Major 

Dinescu et al. ( I2005> have recently re-measured the absolute 
proper motion at the center of the Canis Major over-density, 
which was first determined by Momany et al. ( I2004bl using 
UCAC2 catalogs 9 ). The authors find that while CMa has an in- 
plane rotation (similar to the mean of thick disk stars) it shows 
significant 3cr motion perpendicular to the disk. In turn, they 
find this inconsistent (at a 7cr level) with the expected motion 
of the warp at these Galactic locations (estimated in Drimmel, 
Smart & Lattanzi 2000). This incompatibility lead Dinescu et 
al. ( I2005i to conclude that CMa is part of a satellite galaxy 
remnant. 

To resolve the claimed incompatibility we note that the 
Dinescu et al. ( 2005 1 conclusions rely on two points. The first is 
the assumption that their selected sample [used to measure the 
proper motion and consisting of "likely main sequence stars" 
or blue plume stars] is a photometrically clean sample of gen- 
uine CMa stars, un-contaminated by Galactic disk stars. This 
assumption however is in complete contradiction with the re- 
cent finding of Carraro et al. ( I2005t . who have shown that the 
distribution of the blue plume population in the background of 
30 open clusters (in the Third quadrant and at R G > 7 kpc) fol- 
lows the expected pattern of the Norma-Cygnus (outer) spiral 
arm remarkably well. Thus the CMa over-density resides on an 
outer spiral arm, and the young blue plume population (seen in 
various CMa diagrams) is simply not associated with it. This is 
a particularly important issue because it explains the difference 
in the proper motion at the CMa-center obtained by Dinescu 
et al. ( 120051 using MS stars) and Momany et al. ( I2004bl using 
RGB stars) which in turn brings us back to the situation where 
the proper motion of CMa is indistinguishable from surround- 
ing disk stars 10 . 

9 See however Momany & Zaggia J2005I reporting the presence of 
possible systematic errors in the UCAC catalog. 

10 Note also that in calculating W, Dinescu et al. (i) apply a mean 
V rai ; value referring to the artificially produced V rac i peak at 109 ± 4 
km/s by Martin et al. (2004}; (ii) assume a distance of 8. 1 kpc, differ- 
ent from the most recent determination (7.2 kpc) by Bellazzini et al. 
1 2006 1; and (iii) combine the proper motion based on blue plume stars 
with V md values that are based on red clump stars. 
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Secondly, the fact that the negative W of CMa is irrecon- 
cilable with the expected signature of the Galactic warp (posi- 
tive W) is simply not true, since these expectations are model- 
dependent. We will demonstrate this by using the same cited 
article of Drimmel et al. J2000> . paying more attention to what 
these authors report at the end of their article. 

We first remind the reader that Drimmel et al. ( 120001) used 
Hipparcos OB stars as tracers of the Galactic warp. Therefore, 
it is particularly important to recall recent problems in the in- 
ferred Hipparcos parallaxes of O stars. Schroder et al. (|2004i 
using a sample of 153 stars, conclude that their absolute 
magnitude (calculated from their apparent magnitude and the 
Hipparcos parallax) appear to be much fainter than expected, 
by up to 5 magnitudes. Such huge differences were found to 
arise from the distances at which O stars were located, and 
the level of precision of the parallax measurements achieved 
by Hipparcos. Their Figure 2 clearly shows that large magni- 
tude differences are expected when using the relatively uncer- 
tain Hipparcos parallaxes at distances of already ~ 1 kpc. 

Accounting for this fact, we reconsider the work of 
Drimmel et al. J2000> . In their Figure 11, they report a com- 
parison of the measured vertical velocity component (W) vs 
Galactocentric distances. At Rcc -8.5 kpc, a divergence be- 
tween measured W values of distant OB stars (increasingly 
negative W values) and the expected observable signature of 
the warp (increasingly positive W values) occurs. It remains, 
however, that Drimmel et al. ( 120001 their Section 7) investigate 
the effect of bias, amplitude and precision on the inconsistency 
between the observed and the expected warp W values. In their 
Figure 12, the authors unambiguously show that "negative ver- 
tical motions are finally achieved', if they allow for: (i) a warp 
that has half the amplitude with respect to the one they derive; 
(2) a warp precessing at -25 km s _1 kpc -1 ; and (iii) a 0.5 mag- 
nitude error is added. 

Now, leaving aside the excessively high precision rates 
which Drimmel et al. (120001 also find questionable, it is a mat- 
ter of fact that the remaining two bias strongly contribute in re- 
producing negative vertical velocities. On the one hand, the re- 
sults presented in this paper (see in particular Fig. 2 of Yusifov 
2004 1 show that the warp amplitude as derived by Drimmel et 
al. (FJPOOl is very high. Indeed, at a Galactocentric distance of 
10 kpc, Drimmel et al. note that their warp amplitude is more 
than twice the gas warp amplitude, and not surprisingly they 
halve it. In particular, the authors note that the effect of mis- 
placement of a star to larger distances is that its measured rela- 
tive vertical motion will be smaller than its true relative vertical 
motion. 

On the other hand, there is now enough evidence to suggest 
that the Hipparcos data can have significantly higher errors as 
suggested by Schroder et al. [see also Soderblom et al. (|2005 1 
on a recent confirmation of errors in the Hipparcos parallaxes 
already at the Pleiades distance]. Thus, introducing a 0.5 mag 
(as done by Drimmel et al. 2000) could have been seen as ex- 
treme at that time. Nowadays, however, with confirmed errors 
of more than 2 magnitudes (for the same sample of hot stars) at 
much closer distances, the evidence suggests that the expected 
warp signature can be and is compatible with negative vertical 
velocity. 



Since Dinescu et al. mistake main sequence stars belonging 
to a Galactic spiral arm for genuine CMa population, and given 
the above discussion on the compatibility of negative W values 
with the warp at CMa location, we conclude that the proper mo- 
tion of CMa does not show any peculiar signature with respect 
to Galactic disk stars. 

Appendix D: The "narrowness" of CMa main 
sequence 

Recently Martmez-Delgado et al. (2005) presented a deep 
wide-field B, (B-R) color-magnitude diagram of the CMa cen- 
ter and derived the line-of-sight extent (or depth) of CMa. For 
this task, Martmez-Delgado et al. estimated the observed width 
of the CMa main sequence (<x* „ ,) by selecting stars in a nar- 
row color range, 1 .5 < (B-R) '< 1.55, at B^ 22.1. Their best 
fitting of cr* s tom] (converted into kilo-parsec) yielded FWHM= 
1.95 kpc. Moreover, and for a proper comparison with other 
studied Local Group dwarf spheroidals, they also estimated the 
line-of-sight half-brightness radius (n/2), obtaining ~ 1 kpc. 

Martmez-Delgado et al. (2()()5i noted that this size is "sig- 
nificantly bigger than that of several dwarf galaxies in the 
Local Group", e.g. the n/2 of the Fornax dwarf spheroidal 
is 0.33 kpc (see Table 4 of Irwin & Hatzidimitriou I1995> . 
Nevertheless, one finds that this rather large value found by 
Martmez-Delgado et al. d2005> is often cited in support of a 
"narrow extent line-of-sight" of CMa (Martin et al. 12005 1 Conn 
et al. I2005i or even evidence of a "typical size of a dwarf 
spheroidal galaxy" (Bellazzini et al. l2006t . 

In regards to this, and in light of recent finding of the 
presence of the Norma-Cygnus (outer) spiral arm in the third 
quadrant (Carraro et al. 120051 . we compare the FWHM of 
CMa as obtained by Martmez-Delgado et al. J2005I with typi- 
cal values for Galactic spiral arms. Following Bronfman et al. 
(2000), we estimate the FWHM of the most distant spiral arm 
at R/R = 1.3, where R =8.5 kpc. Relying on their Fig. 8, we 
estimate a FWHM of ~ 1 .7 kpc, very much in accordance with 
the Martmez-Delgado et al. ( I2005i 1.95 kpc value for CMa. 
Indeed, one must also bear in mind that an outer spiral arm 
would also be more flared, thereby a FWHM of ~ 1.7 kpc sets 
only a lower limit. One last piece of evidence that the reported 
CMa FWHM is compatible with a spiral arm feature comes 
from the study by McClure-Griffiths et al. (120051 . The detected 
H 1 arm in the Fourth and Third quadrants appeared to be ap- 
proximately 1 to 2 kpc thick along the line of sight. 

We therefore conclude that the CMa depth as measured by 
Martmez-Delgado et al. ( I2005i is more compatible with a dis- 
tant spiral arm rather than a dwarf spheroidal. 

Appendix E: CMa and the Orion arm connection 

Recently, Moitinho et al. (120061 used young (~ 100 Myr) 
open clusters and field "blue plume" population to reconstruct 
the spiral structure in the Third Galactic quadrant. They con- 
firmed the detection of a structure interpreted as a natural 
extension of the Norma-Cygnus spiral arm. Their conclusion 
(based on U BVRI photometry of OB objects) is consistent with 
our finding of an extended MW disk (using 2MASS M-giants). 
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Moitinho et al. went a step further and proposed that the CMa 
M-giant over-density is simply the result of looking along the 
extension of the local Orion arm in the Third quadrant. 

In regards to this, we note that the Moitinho et al. {2006 1 
and Carraro et al. (2005? papers have been valuable in clear- 
ing the connection between the blue plume population (erro- 
neously attributed to CMa by Martmez-Delgado et al. 2005 
Bellazzini et al. 1200 4 and Dinescu et al. 12005 1 and both the red 
clump and red giant populations (that first brought to the detec- 
tion of the CMa over-density). Nevertheless, the Moitinho et 
al. analysis does not address, nor preclude, the very occurrence 

of the CMa over-density at b 8°. Indeed, the detection of a 

spiral arm composed by young blue objects (best identified in 
U BV photometry) is not sufficient to affirm/conclude a Galactic 
nature for the un-ambiguous presence of the older (2-10 Gyr) 
CMa M-giant over-density (best seen in the 2MASS infrared 
photometry). The asymmetry in the distribution of the MW 
disk stars (Galactic warp) remains the explanation for the ob- 
served asymmetric number density at b ~ -8° and b ~ +8°. 
The warp signature is evident in Fig. 2 of Moitinho et al. 
(120061 



